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1 •  INTRODUCTION 


Gun  accuracy  essentially  encompasses  the  complete  delivery  mechanism  of 
one  or  many  projectiles  fired  from  a  gun  to  a  target  intercept  space-time 
volume  so  that  the  warhead  will  be  able  to  defeat  the  target.  The  main  thrust 
of  the  ongoing  US  work  is  directed  toward  tank  gun  accuracy  and  projectile 
launch  integrity  in  order  to  increase  the  hit  probability  of  our  tank 
fleets.  Much  of  this  involves  product  improvement  of  fielded  weapons  and 
development  support  for  new  systems.  Hence,  most  of  the  work  involves  system 
engineering  and  analysis  and  is  hardware  related  and,  as  such,  is 
classified.  The  Ballistic  Research  Laboratory  (BRL)  is  in  the  process  of 
actively  exploiting  innovative  ideas  and  concepts  to  assure  that  projectiles 
will  leave  the  muzzle  in  the  proper  aim  direction  by  reducing  the  gun  related 
variation  of  projectile  launch  parameters  to  that  obtained  when  firing  from 
precision  mounts.  To  support  these  tasks  and  develop  the  appropriate 
technology  base,  BRL  is  pursuing  a  basic  research  program  which  addresses, 
among  other  technology  areas,  the  projectile  launch  from  the  gun. 

In  particular,  we  are  concerned  with  the  sources  and  mechanisms  which 
cause  the  bias  of  the  impact  point  relative  to  the  target  as  well  as  with  the 
errors  and  impact  distribution  due  to  parameter  variations.  For  all  practical 
purposes,  the  effect  of  exterior  ballistic  related  sources  as  well  as  the 
propagation  of  variations  in  the  launch  parameters  on  the  projectile  flight 
trajectory  are  well  understood  and  predictable  with  sufficient  accuracy. 
Contrarily,  the  sources  and  processes  which  determine  the  variation  of  the 
projectile  launch  parameters  are  least  understood.  Though  a  reasonable 
qualitative  understanding  exists,  quantitative  relationships  remain  a 
problem.  Unfortunately,  a  methodology  correlating  projectile  and  gun 
parameters  with  the  projectile  dynamics  immediately  after  projectile  release 
from  the  muzzle  and  satisfying  the  needs  of  quantitative  analysis  has  not  been 
developed  yet.  To  remove  this  obstacle  and  to  develop  the  required 
methodology  for  a  quantitative  simulation  of  projectile  launch,  BRL  is 
pursuing  a  multi-pronged  approach  including  both  theory  and  experiment.  The 
experimental  program  addresses  both  the  validation  of  gun  system  accuracy 
models  under  development  and  the  improvement  of  our  diagnostic  capability. 

One  part  focuses  on  the  investigation  of  the  disengagement  mechanism  of  the 
projectile  from  the  cannon.  This  includes  the  adaptation  of  available 
instrumentation  and  measurement  techniques  along  with  the  development  of  new 
ones  suitable  for  the  measurement  of  the  muzzle  motion  with  the  appropriate 
temporal  and  spatial  resolution  and  without  any  interference  to  the  projectile 
and  gun  dynamics.  The  muzzle  dynamics  is  intrinsically  connected  to  the 
projectile  motion,  i.e.,  muzzle  aim,  transverse  velocity  and  angular  rotation 
during  shot  exit  are  manifested  in  and  reflect  the  projectile  launch 
parameters . 

Most  of  the  instrumentation  and  measurement  techniques  for  monitoring  gun 
barrel  motion  are  directly  applicable  to  the  measurement  of  the  six  degrees- 
of-freedom  motion  of  the  muzzle  section  of  the  gun  barrel.  A  brief  review  of 
barrel  motion  measurement  techniques  is  provided  in  Appendix  A. 

Unfortunately,  there  are  certain  difficulties  involved  in  selecting  the 
appropriate  measurement  techniques  and  combining  them  into  an  instrumentation 
setup  which  allows  the  extraction  of  the  six  degrees-of-f reedom  motion  of  the 
muzzle  with  the  required  accuracy  from  the  observables.  For  instance,  the 
effect  of  longitudinal  and  radial  stress  waves  produced  in  the  tube  by  the 
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traveling  projectile  and  combustion  gas  pressure  as  well  as  the  abrupt 
decorking  of  the  tube  is  largest  at  the  muzzle  end  and  manifests  itself  in  the 
recoil  motion  of  the  muzzle  as  the  projectile  moves  through  it.  We  have 
observed  short  time  reductions  of  the  muzzle  recoil  velocity  by  more  than  one 
fourth  from  its  value  prior  to  shot  exit.  Therefore,  the  recoil  motion  must 
be  measured  directly  at  the  muzzle  and  cannot  be  inferred  from  measurements  at 
other  locations.  Also,  because  of  the  high  frequency  content  of  these  stress 
waves,  the  use  of  piezoelectric  accelerometers  mounted  on  the  outer  surface  of 
the  cannon  may  not  be  expedient  for  the  measurement  of  the  axial  and  the 
transverse  muzzle  motion.  Other  instrumentation,  such  as  the  electrooptical 
or  electromagnetic  displacement  sensors,  must  be  employed  which  do  not 
interfere  with  the  projectile  or  the  tube  motion  and  are  not  affected  by  the 
radial  vibrations  of  the  tube  and  the  precursor-blast  environment.  In 
addition,  a  time  resolution  on  the  order  of  microseconds  or  less  and  a  spatial 
resolution  on  the  order  of  microns  are  desirable  to  adequately  describe 
projectile  exit  phenomenology. 

The  Schmidt  displacement  transducer  ( SDT)  instrumentation  and  measurement 
technique  recently  developed  at  BRL  for  the  determination  of  the  projectile 
motion  during  shot  exit  [1] [2] [3] [4]  can  also  be  applied  to  the  measurement  of 
the  muzzle  motion.  This  displacement  transducer  can  provide  data  which,  after 
the  appropriate  data  analysis,  allow  the  determination  of  five  parameters  of 
the  six  degrees-of-freedom  motion  of  the  muzzle.  Only  the  rotation  about  the 
bore  axis  cannot  be  determined.  The  technique  has  been  applied  and  checked  in 
a  recent  experiment  [5]  [6]  and,  subsequently,  has  been  employed  in  a  tank  gun 
accuracy  firing  test  to  obtain  the  muzzle  aim  at  shot  exit. 

This  report  contains  a  brief  description  of  the  SDT  (Section  2.)  and 
discusses  the  data  analysis  concept  and  procedure  (Section  3.)  which  form  the 
basis  of  the  computer  program  (Section  4.).  The  data  used  for  illustrating 
the  analysis  program  are  from  the  cited  experiments. 

2.  SDT  INSTRUMENTATION  AND  MEASUREMENT  TECHNIQUE 

The  Schmidt  displacement  transducer  for  monitoring  transverse 
displacements  (Figure  2.1)  as  well  as  its  function  are  described  in  Reference 
[4].  As  shown,  the  SDT  consists  of  two  pairs  of  semicircular  inductors 
peripherally  connected  to  an  oscillator  and  concentrically  arranged  in  a  ring 
configuration,  ninety  degrees  apart.  This  ring  is  placed  concentric  with  and 
perpendicular  to  the  tube  axis  at  the  desired  measurement  location,  thus 
establishing  an  electromagnetic  coupling  between  the  active  sensor  and  the 
conducting  tube  surface.  The  currents  from  each  sensor  pair  are  fed  to  a  low- 
gain  differential  amplifier  and  its  output  is  then  amplified  to  a  suitable 
level  for  recording.  The  polarity  and  magnitude  of  these  currents  are  a 
direct  measure  of  the  direction  and  the  magnitude  of  the  closeness  of  the 
outer  barrel  surface  to  the  two  semicircular  inductors  and,  when  calibrated, 
yield  the  appropriate  transverse  displacement  of  the  tube  relative  to  the 
electromagnetic  center  of  the  SDT  ring.  Employing  two  or  more  SDT  systems  at 
the  muzzle  allows  the  pointing  vector  of  the  muzzle  to  be  extracted  from  the 
displacements.  The  functional  relationship  between  the  measured  voltages  and 
the  corresponding  transverse  displacements  has  been  found  to  be: 
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Figure  2*1.  Sensor  Arrangement  And  Block  Diagram  Of  Electronic  Circuit  Of  The 
SDT 


V 

x 

x  =  - -  ■  1  1  " 

a  yl  +  aV2  +  aV2 
o  xx  y  y 

where  Vx>  and  V  denote  voltage  and  x  and  y  the  corresponding  displacements. 

The  coefficients  occurring  in  these  relations  will  vary  from  setup  to  setup 
and  must  be  determined  prior  to  each  firing  experiment.  Since  the  sensor 
rings  are  considerably  lighter  than  the  cannon,  the  calibration  is  best  done 
by  mounting  the  sensor  ring  on  positioning  devices  and  displacing  the  sensor 
with  respect  to  the  resting  cannon. 

The  axial  or  z-cornponent  of  the  tube  motion  can  also  be  obtained  by 
employing  the  SDT  technique.  This,  however,  requires  that  the  exterior  cannon 
surface  be  modified.  This  is  achieved  by  alternating  spatially  either  the 
radius  or  the  electric  conductivity  of  the  outer  surface  of  the  cannon. 

Though,  in  principle,  the  sum  signal  from  the  SDT  system  can  be  used  directly 
for  monitoring  the  translation  in  the  z-direction,  it  is  advisable  to  separate 
the  measurement  of  the  z-displacement  from  that  of  x  and  y.  Two  sensor 
designs  in  which  unwanted  contributions  from  the  transverse  and  rotational 
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tube  motion  are  sufficiently  suppressed  are  practicable:  one  consisting  of 
four  individual  sensor  loops  ninety  degrees  apart,  aligned  parallel  with  the 
tube  axis  with  the  four  signals  summed  (Figure  2.2  (a))  and  the  other 
consisting  of  a  single  loop  arranged  concentrically  around  the  tube  (Figure 
2.2  (b)).  The  relationship  between  the  recorded  voltage  and  the  corresponding 
displacement  is  given  by 
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Figure  2.2.  Sensor  Arrangement  And  Signal-Displacement  Correlation  For  The 
SDT 
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where  is  the  accumulated  recoil  displacement  up  to  the  extremum  E^+j  and 

AZ,  the  distance  between  the  extrema  and  E^.  The  coefficient  of  the 

polynomial  P  (C)  may  vary  from  interval  to  interval  and  must  be  determined 
prior  to  the  firing  experiment.  The  width  of  the  periodic  bands,  w,  depends 
on  the  design  and  the  accuracy  requirement.  For  the  four  sensor  loop  design  w 
should  be  a  full  or  one-third  sensor  length  and  for  the  single  sensor  loop 
design  w  should  be  about  the  mean  distance  from  the  sensor  ring  to  the  cannon 

surface.  From  the  analysis  point  of  view,  the  highest  accuracy  and  linearity 

should  occur  during  shot  exit.  This  requires  that  the  initial  location  of  the 
sensor  at  shot  start  be  selected  such  that  its  location  at  shot  exit  is  in  the 

vicinity  of  sin  P(0  «  0  or,  equivalently,  the  band  edge  should  be  close  to 

the  center  plane  of  the  sensor.  Since  radial  tube  expansion  and  contraction 
may  slightly  shift  the  signal  level,  w  should  be  short  enough  to  include  the 
recording  of  the  extrema  enclosing  the  measurement  interval  of  interest,  thus 
allowing  self-adjusting  conversion  of  the  signal. 

Special  care  has  to  be  taken  in  the  mounting  arrangement  for  the  SDT 
rings  in  order  that  neither  the  muzzle  blast  nor  the  counterrecoiling  muzzle 
will  damage  the  instrumentation.  In  firing  tests  with  tank  guns  it  has  been 
found  that  it  is  best  to  physically  disconnect  the  sensor  rings  from  the 
mounts  after  the  projectile  has  left  the  muzzle  and  let  them  freely  move  with 
the  recoiling  gun  '.  For  small  caliber  guns,  this  provision  is  not  necessary. 

3.  DATA  ANALYSIS  GONCEPT  AND  PROCEDURE 

As  pointed  out  earlier,  the  desire  for  gun  accuracy  is  the  driving  force 
behind  most  of  our  cannon  motion  measurements.  These  measurements  must 
provide  data  which  allow  the  correlation  of  target  location,  aiming  point  of 
muzzle  after  laying  the  gun  (shot  start),  aiming  point  and  transverse  as  well 
as  angular  velocity  of  muzzle  at  shot  exit,  initial  free  flight  motion  of 
projectile,  and  target  impact  location.  The  obtaining  of  appropriate  muzzle 
motion-related  parameters  for  direct  fire  weapons  mounted  on  resting  platforms 
is  addressed  here. 

3*1.  Experimental  Arrangement 

It  is  presumed  that  the  transverse  displacements  of  the  muzzle  section  of 
the  cannon  is  monitored  at  m  locations,  where  m  is  at  least  two,  and  the 
muzzle  recoil  at  mf  locations,  where  mf  is  at  least  one,  with  sensors  which 
measure  these  displacements  with  respect  to  the  ground  (Figure  3.1).  An  m=3 
and  mr  =  2  system  is  currently  being  developed  at  RRL,  consisting  of  three 
transverse  SDT  units  about  three  bore  diameters  apart  and  two  axial  single 
loop  SDT  units  located  in  between  to  measure  five  of  the  six  degrees-of- 
freedom.  Only  the  rotation  of  the  tube  about  its  main  axis  is  not 
obtainable.  Though  the  data  analysis  is  general,  it  has  been  developed  with 
this  particular  case  in  mind. 


*  ) 

7  An  upcoming  BRL-MR  by  J.Q.  Schmidt  and  T.L.  Brosseau  will 
describe  a  mounting  arrangement  for  the  SDT  rings  which  was 
successfully  employed  in  a  recent  tank  gun  accuracy 
experiment . 
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m  =  3  m  -  2  m=l 


2(  =  Y) 


Figure  3.1.  Sensor  Locations  For  Monitoring  Muzzle  Motion 


It  is  further  assumed  that  the  instrumentation  is  completely  decoupled 
from  the  gun  and  its  platform;  i.e.  ,  work  on  the  gun  can  be  done  without  any 
disturbance  of  the  sensors.  The  frame  of  reference  provided  by  the  array  of 
sensors  must  stay  the  same  from  the  time  the  gun  is  layed  until  the  shot  has 
left  the  muzzle.  This  requirement  allows  one  to  relate  the  position  of  the 
tube  axis  at  the  individual  measurement  locations  to  a  reference  line 
established  prior  to  the  actual  firing  and  to  calibrate  the  muzzle-pointing 
vector  as  obtained  from  the  sensors  to  the  aiming  point  established  by  bore 
sighting  (Appendix  B) . 

3.2.  Analysis  Procedure 

It  is  presumed  that  the  recorded  signals  have  been  converted  into 
voltages  by  accounting  for  the  appropriate  calibration  for  the  total 
instrumentation  and  recording  chain  and  are  arranged  into  time  series  with  a 
common  time  basis, 


V( v ,  y,  X)  ,  v  =  1,2  ,  p  =  1,2,  ...  m  ,  X  =  1,2, 

V( 3 ,  p'  ,  X),  p'  =  1,2,  ...m'  ,  X  =  1,2, 

where  X  is  the  index  for  the  discrete  time  t  =  X  At , 
and  H  is  the  length  of  the  available  data  window. 

3.2.1.  Data  Conversion 

These  voltages  are  then  converted  into  displacements  using  the 
appropriate  mapping  correlation.  For  the  transverse  SDT  signals  we  have, 
according  to  Eq.  (2.1), 


A  , 
A  > 


(3.1) 


At  is  the  time  increment , 
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S(v,  p  ,  A)  =  V(v,  p,  A)/[  C. 


V7T 

v  o 


CvlV<1- 


X)  +  Cv^<2,  u,  X)z  8 


*) 

S(v  ,  p  ,  X)  =  S(v  ,  p  ,  A)  -  Sq(v  ,  p)  , 

v  =  1,2,  p  =  1,2,...  m ,  A  =  1,2,  ...  £ 


(3.2) 


where  S  (v ,  p)  is  an  alignment  function  as  defined  in  Appendix  B.  The 
locations  along  the  tube,  Z(p)  =  Z(p ,  X),  where  these  transverse  displacements 
are  recorded,  move  with  the  tube  axially  and  are  computable  from  the  axial  SDT 
signals.  Their  conversion  into  displacement,  though  more  complex,  is  also 
straightforward  if  one  follows  the  procedure  sketched  below. 


V  p’,  p’  e  [1,  2, 


0  find  the  extrema  of  the  time  series  V(3,  p',  X),  X  =  1,2,... ft,  and  their 
corresponding  times 


E(p',  j),  ft(j)  =  ,  j  =  1,2, 


0  ¥  j,  j  e  [1,  2,  ...J^,]  : 


0  for  j  =  1  : 

identify  the  z-interval  to  which  E(p',l)  belongs  i  =  i  , 
let  AE  =  E.-E.  ,  where  E.-E.  is  taken  from  calibration,  and 

a  1  1  i  1  X  It  i 

set  A  =  E(p',l)  -  A  E,  Zq  =  0,  A  =  0 

0  let  AZ  =  Z  -  Z  where  Z.-Z. , .  is  given  by  the  calibration,  and 

1  1x1 f  X  111  A 

set  B  =  A,  A  =  E(p ' ,  j)  ,  Aq  =  A  +  1,  A  =  £(j) 

0  ¥  X  ,  X  e[AQ  ,  A]  : 

let  £  =  2[V(3,p  ’,  X)  -  B]  /  [  A  -  B]  -  1 
and  compute 

S(3,  p',X)  =  Z  +AZ[  Vo  +  —  arcsin  P(0] 

O  A  7T 


*)  FORTRAN  type  substitution  statements  are  used  throughout. 

**)  Special  notation  is  being  used  where  the  symbols  have  the  following 
meaning:  ¥...  for  all  and  e...  from  the  interval 
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0 


if  j=l:  subtract  S(3,  y ' ,  1)  from  the  calculated  displacements  to 
equate  the  z-di splacement  at  shot  start  to  zero 

VA,  Ae  [  1 ,  X<  1) 3  :  S(3,y',A)  =  S(3,y',A)  -  S(3,y\l) 

A 

set  ZQ  -  S(3,  p  ,  £(j))  and  x= i  1  to  provide  the  initial  displacement 
and  identification  for  the  next  z-interval  (3.3) 


3.2.2.  Modeling  Of  Muzzle  Flexure 

To  develop  the  above  mentioned  data  analysis  algorithm  which  allows 
the  extraction  of  five  of  the  six  degrees-of-f reedom  mot  ion  parameters  from 
the  time  series,  we  have  to  make  certain  model  assumptions. 

The  tube  flexure  near  the  muzzle  end  (z=0)  may  be  modelled  as  a 
polynomial  of  k-th  order  in  z: 


S(v,  z,  X)  -  £  C(k,  v3  X)zK,  v=1,2,  X*l,2...£  .  (3.4) 

0 

For  the  end  conditions  for  a  free  beam  at  z^O  we  may  assume  zero  bending 
moment  as  well  as  zero  shearing  force.  Thus, 


[EIS11] 


z=0 


0  S"(v5  0,  X)  =  0,  and 


(3.5) 


[EIS" ] 


z=0 


0 


S,M(v,  0,  X)  “  0 


(3.6) 


where  Eq.  (3.5)  is  the  expression  for  bending  moment,  Eq.  (3.6)  the  expression 
for  shearing  force,  E  Young's  modulus  of  elasticity  and  I  =  I(z)  the  moment  of 
inertia.  From  that  we  have  the  requirement  that 


C(2,  v,  A)  -  C(3,  v,  A)  =  0  and  k  <  y+1 


(3.7) 


As  the  gun  barrel  recoils,  the  sensor  locations  with  respect  to  the  muzzle 
face  will  vary  with  the  recoil  as 


z(p,A)  =  zq(m)  +  /  dA  z(z,A)  ,  n  =  1,2, ...m  , 


(3.8) 


where  z(z,A)  is  the  recoil  velocity  of  the  muzzle  which  is  obtainable  from  the 
S(3,  y  ,  a),  y  - l,2,...m'.  Because  these  recoil  displacements  are  not 


recorded  at  the  locations  z(y),  y=l,2,...m  ,  and  because  the  tube  response  to 
the  travelling  projectile  and  combustion  gas  pressure  is  large  and  transient 
at  the  muzzle,  straightforward  linear  inter-  or  extrapolation  may  not  be 
expedient,  especially  if  the  muzzle  motion  is  to  be  determined  with  very  high 
accuracy. 

The  recoil  velocity  at  the  locations  z(p')  may  be  computed  by  convoluting 
the  time  series  S(3,  p 1 ,  A)  with  an  appropriate  FIR  (finite  impulse  response) 
differentiator  D(x')  [7], 


S(y',A)  =  E  D(X’)S(3 ,uf,  A“A 1 )  ,  p,=l,2,...  m'  ,  X=l,2...Jt  .  (3.9) 

A  1 

This  velocity  history  (Figure  3.2)  consists  of  a  trend  describing  the  gross 


OBSERVED  RECOIL  VELOCITY 


* 

ZT 


Zsw 

A 


SHOT  EXIT 
X"  \+ 


SQUARE  WELL-LIKE  MODULATION 


Figure  3.2.  Muzzle  Recoil  Velocity 

recoil  motion  of  the  tube  at  that  particular  location  and  a  square  well-like 
modulation  derived  from  the  transient  in-bore  loads.  The  arrows  show  the 
behavior  of  the  well  as  the  measurement  location  is  moved  towards  the  muzzle 
face:  the  width  enveloping  the  shot  exit  time  narrows  and  the  depth 

increases.  The  left  wall  correlates  to  the  passage  of  the  rotating/driving 
band  of  the  projectile  travelling  towards  the  muzzle  and  the  right  wall  to  a 
tube  stress  relaxation  wave  produced  by  the  decorking  of  the  pressurized  tube 
and  travelling  from  the  muzzle  back  to  the  breech.  To  separate  the  gross 
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recoil  velocity  and  the  well-like  modulation  in  the  presence  of  noise,  we  may 
proceed  in  the  following  way: 


0  V  p 1 ,  p 1  e  [ 1,2,  . .  .m'  ]  : 

0  Compute  an  upper  and  lower  envelope  of  the  data,  their  difference, 
and  their  maximal  and  mean  values: 


A  S  -  A  &  —  0 

max 


V  A  ,  Xe  [AA  +  1 , . . .£-AA ] ,  where  2  AX  +  1  is  the  arc  length. 


&+(p 1 5X )  =  max{ 5(p ' ,A“AA) , • • . S(y 1 ,A+AA )}  .. .upper  bound  , 
&  (y',A)  =  minj  §(y  * , A~AA ) , • • • §(y 1 ,A+AA ) }  ...lower  bound  , 


AS(pf,A)  -  S  (  y  f  ,  X  )  ^  S  (  y  f  » X  ) 


. . .difference  , 


A&  =  maxfA§  ,  A&Cu^A)}  ...maximum  ,  and 

max  1  ttiav  p  i 


max 


A&  -  A&  +  A§(y  1  ,A)  , 


AS  =  AS/U-2AA) 


• • .mean  • 


The  behavior  of  the  resulting  function  is  illustrated  in  Figure  3.3. 
S  AS 
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0 


.  .  - — f-  +-  +  + 

Determine  the  times  A  ,A  ,A  ,  and  X  at  which  the  difference 

curve  intersects  the  line  =  (AS  +  AS)/2  and  compute  from  them 

max 

the  time  locations  of  the  left  and  right  well  boundaries, 
l  (p1)  =  V2  (X  +A  +)  and  t+(y')  «  V2  U+~+A++)  . 

0  Using  spline  fitting  or  least  squares  (LSQ)  technique,  set  forth  a 
function  (e.g.,  a  sequence  of  cubic  splines  or  a  higher  order 
polynomial)  such  that 

f(p'3X)  L=^  {SCp'jA)}  for  X  <  £  (y 1 )  and  X  >  £+(p')  , 

where  f(y',X)  now  describes  the  gross  recoil  velocity.  Subtraction  of  this 
function  from  the  data  yields  a  time  series  which  contains  the  leftover  noise 
and  the  well-like  modulation: 

AS(y 1 ,£ )  =  S(p',A)  ~  f(p’,X)  ,  X  =  1,2, . . .£. 

0  Using  a  LSQ  procedure,  the  depth  of  the  well  can  be  obtained 

g(u',X)=def  g0+glX  L"Q  {AS(p  *  ,  A  ) }  ,  Xen~(u')+6fc,Jl  +  (u')-6£]. 

0  With  that,  the  decomposition  of  the  individual  velocity  signal  into  a 
gross  recoil  velocity  and  a  time  dependent  well-like  modulation  is 
completed , 


§(u',x)  =  f(p  '  ,A)+e(x-*  (j!  '))e(ji+(y)-x)g(y  '  ,\)  +  n(u',x)  , 


where  e(x)  0  for  x^O  and  e(x)  —  1  for  x>0  and  n(p 1 ,x)  is  the  remainder  after 
subtracting  the  gross  recoil  velocity  and  the  well  from  the  data. 

0  From  these  individual  functions  a  new  function  describing  the  recoil 

motion  of  the  muzzle  relative  to  the  instantaneous  measurement  location 
z(y,=l,X)  can  be  derived: 


z(?3X)  -  F(^,x)+e(x~L  (^))e(L  ( £ )~A )G( £ ,A ) ,  where 

k' 

F(C,i)  -  de£  I  L-Q  {«„',»>)  , 
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!r  » 


G(4  ,X) 

"  def 

2  g 
0 

k' 

L_(4)  = 

def 

Z 

o  K 

k' 

L+(4)  = 

def 

Z 

~  K 

.K  LSO 


s  V*  L|°  { £~(u  ')}  , 


Z  LIQ  {  £+(u  •  )}  , 


for  k'CU'-l),  c  =  z-z(y'  =  l,X)  and  4  (u  '  )=ZQ(y  '  )-z  ( 1)  .  (3.10) 


By  continuously  updating  the  location  of  the  muzzle  we  can  obtain  the 
translation  of  this  location  relative  to  the  resting  frame  of  reference: 

^  z  Zo ^ ^  ^  —  z(£,l)  ...  initial  location  and  velocity  with  respect 


to  2  (y  1 
0 

X ,  Xe [2, . .  .£]  : 

4  =  4+&X4 

• • .updated 

?  =  z(4,X) 

• . .updated 

z  (X)  =  z  (1)  +  4 

,  .  0 

• . .updated 

z(X )  =4 

•  •  .updated 

Setting  z  =  ZQ(y),  y=l,2,...m,  provides  the  cannon  recoil  motion  at  the  sensor 
locations  which  monitor  the  transverse  muzzle  displacements.  Insertion  of  the 
appropriate  quantities  into  Fq.  (3.4)  yields  a  set  of  algebraic  equations, 


V  v  and  X  ,  v=l , 2  and  X  e[l, 2, 


1  z(l,X)  z4(l,X)...zk(l,X) 

C(0,v,X) 

S(v,l,X) 

1  z(2,X)  z4(2,X)...zk(2,X) 

C(1,V,X) 

S(v  ,2 ,X  ) 

•  *  . 

» 

C(4,v,X) 

= 

*  .  « 

1  z(m,A)  z  (m, A  )  . . . (m,A  ) 

C(k,v,X) 

S(v  ) 

(3.12) 


from  which  the  coefficients  C(k,v,X),  k=0,  1 ,4, . .  .k«y  +  l) ,  can  be  computed 
either  directly,  if  k=y+l,  or  via  LSO,  if  k<y+l. 
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3. 2 .3.  Translation  And  Rotation  Parameters 

The  translation  of  a  point  lying  on  the  muzzle  axis  with  respect 
to  the  ground  frame  of  reference  is 

k 

S(v,z,A)  =  Z  C(k,v,\)  z(A)K  ,  v=l  ,2  , 
o 

S(3,z,A)  =  2(A)  ,  A  =  1,2,.. •£  .  (3.13) 


One  can  now  obtain  its  linear  transverse  velocity  by  convoluting  the  time 
series  with  an  appropriate  FIR  differentiator: 

S(v,z,A)=  ^S(v,z,A)  =  E  D(A')S(v,z,A-A’),  v-1,2,  A=l,2,...£. 

A  * 

(3.14) 

The  linear  axial  velocity  S(3,z,A)  =z(z,A)  is  given  by  Eq.  (3.11). 

The  slope  of  the  muzzle  at  position  z  is  given  by 

d  k  K  i 

Sf(v,z,A)  =  jz  S(v,z,A)  =  Z  k  C(k,v,A)z(A)  ,  v-1,2,  X~1,2,...A. 

o 

(3.15) 

From  Eq.  (3.15)  and  neglecting  any  rotation  of  the  tube  axis  we  can  determine 
the  axes  of  the  moving  coordinate  system  which  is  located  at  the  axial 
location  z  (Figure  3.4)  relative  to  the  ground  frame  of  reference 


eJ(A)  «  I  (1,0,  -  sj)  , 

^2(X)  =  b  (-SiS2’  1+s'l2’  ~sp  ’ 

epA)  =  i  (sj,  sp  1)  ,  where 


* 

u  =• 

Vl+sj2  ,  v  =  Vl  +  sj2  +  s^2  , 

sr 

S' ( 1  ,z ,A  ) ,  s£  =  S' (2,z,A )  ,  and 

A  = 

1,  2,  .. .  1  . 

(3.16) 

The  axis  given  by  e^  may  also  be  described  by  two  successive  rotations 
6  and  <j)  relative  to  the  reference  axis.  The  expression  for  0  and  <}>  in  terms 
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of  the  instantaneous  tube  axes  are 


0(A)  =  arcsin  /\/l+e^  +  e32  1  » 


4>  ( A )  =  arccos  [  e32^e31  +  e32  J  * 

where  e^  =  e^(A)  >  e32  =  e32^X)  ,  A=l,2,...£  .  (3.17) 


3'  (INSTANTANEOUS 
TUBE  AXIS) 


3  (REFERENCE  AXIS) 


Figure  3*4*  Rotation  Of  Tube  Axis  With  Respect  To  Ground  Reference  System 

The  rate  of  change  of  0  and  <j>  may  be  determined  by  numerical  differentiation: 


0(A)  =  Z  D(A')6(A-A')  and  $(a)  =  Z  D(  A  ' )  <f»  ( A~A ' )  .  (3.18) 

A '  A' 

Having  determined  the  translation  and  rotation  of  a  coordiante  frame 
centered  on  the  muzzle  axis  we  can  express  any  vector  given  with  respect  to  a 
moving  coordinate  system  as  a  vector  with  respect  to  the  ground  frame  of 
reference  as 


3 

V(v,X)=S(v  ,z,x)  +  £  V*(z;vf  ,X)e' t  ,  v=sl,2>3,  X  =  1 , 2 , .  •  .£  . 

1  v  v 

(3.19) 

It  should  be  noted  that  the  derived  quantities  depend  on  the  initial  gun 
curvature  which  is  very  seldom  obtained  experimentally.  Equating  S  (v>y), 
v=l,2  and  y=l,2,...m,  in  Eq.  (3.2)  to  the  mean  value  of  the  displacement  time 
series  before  any  motion  occurs,  the  above  equations  describe  the  temporal 
change  in  the  displacements  with  respect  to  the  unknown  initial  tube  flexure. 
Because  of  the  free  end  conditions  of  the  tube  at  the  muzzle  face  the 


equations  will  yield  the  correct  motion  at  the  muzzle  end  regardless  of  the 
downtube  curvature  before  shot  start. 

4.  OUTLINE  AND  DESCRIPTION  OF  COMPUTER  PROGRAMS 
4.1.  General  Description 

The  BRL  CYBER  computer  system  consists  of  three  mainframes:  mainframe 
A(MFA),  which  is  a  CYBER  750,  mainly  used  for  input/output,  mainframe  B(MFB) , 
a  CYBER  825,  for  computer  graphics,  and  mainframe  Z(MFZ) ,  a  CYBER  7600,  for 
floating  point  computations.  In  addition,  MFA  and  MFB  share  a  common  file 
space  and  both  use  the  NOS  operating  system.  MFZ,  though,  uses  the  SCOPE 
operating  system.  The  commercial  package  DISSPLA  [8]  was  chosen  to  generate 
the  plots.  Version  8.2  of  DISSPLA  is  resident  on  all  mainframes,  while 
Version  9.0  is  available  only  on  MFA  and  MFB.  Plots  are  available  on  a 
CALCOMP  plotter  or  a  printer  attached  to  an  interactive  terminal  (in  this 
case,  a  TEKTRONIX  4695  printer  interfaced  to  a  Tektronix  4107  terminal).  One 
computer  program  uses  the  International  Mathematical  and  Statistical  Library 
(IMSL)  [9],  which  is  also  available  on  all  mainframes. 

The  data  analysis  computer  package  consists  of  two  programs,  MUZM040  and 
MUZPRED,  both  written  in  FORTRAN  IV.  The  first  program  is  a  preprocessor  and 
the  second  one  does  the  actual  computations.  Program  MUZM040  inputs  the  raw 
data  (in  counts),  converts  the  data  into  engineering  units  and  then  filters 
and/or  differentiates  the  data.  The  input  data  are  assumed  to  be  on  a  data 
file.  In  our  case,  a  Nicolet  oscilloscope  is  interfaced  with  a  Hewlett- 
Packard  9845 C  microcomputer  which,  in  turn,  communicates  with  the  BRL  CYBER 
computer  system.  The  data  are  retrieved  from  the  Nicolet  disk,  stored  on  the 
HP9845C  disk,  and  then  transferred  to  a  file  on  MFA.  This  transfer  is 
described  in  Reference  [10].  The  output  data  file  from  the  first  program  is 
then  used  as  an  input  to  the  second  program. 

Listings  of  the  programs,  the  job  control  language  (JCL)  to  run  the 
programs,  descriptions  and  formats  for  the  input  data  and  a  sample  case  along 
with  sample  input  and  output  are  given  in  Appendices  C  and  D. 


Figure  4.1.  Flowchart  For  Program  MUZM040 
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4.2.  Computer  Program  MUZM040 


The  flowchart  for  MUZM040  is  shown  in  Figure  4.1*  The  main  program  calls 
the  subroutines  in  proper  sequence  as  determined  by  the  input  control 
variables:  IPLT1,  IPLT2,  IPLT3,  and  IPRT.  If  the  control  variable  equals  1, 

the  operation  is  done;  otherwise,  not.  IPLT1  controls  the  plotting  of  the 
displacements;  IPLT2  controls  the  plotting  of  the  velocities;  IPLT3  controls 
the  plotting  of  the  accelerations;  and  IPRT  controls  the  saving  of  the 
processed  data  for  future  use. 


4.2.1.  Subroutine  C0EF1 


This  subroutine  inputs  the  coefficient  array  and  other  necessary 
parameters  of  a  previously  generated  FIR  lowpass  or  highpass  filter. 


4.2.2. 


Subroutine  C0EF2 


This  subroutine  inputs  the  coefficient  array  and  other  necessary 
parameters  of  a  previously  generated  FIR  differentiator. 


4.2.3. 


Subroutine  CONVOL 


This  subroutine  does  a  linear  convolution  of  the  NFILT 
coefficients  with  a  data  array  of  length  NFILT.  This  process  then  proceeds  in 
moving  arc  fashion  along  the  entire  data  array.  The  convolutions  at  the 
endpoints  are  accounted  for  by  reflection  of  the  data  at  the  endpoints.  See 
the  listing  in  Appendix  C  for  more  details. 


4.2.4.  Subroutine  CONVT 


This  subroutine  converts  the  displacement  data  from  counts  to 
voltages  and  then  to  engineering  units  (microns,  in  this  case).  The 
conversion  is  nonlinear  as  discussed  previously  in  Section  2  and  as  expressed 
in  Eq.  (2.1).  The  data  are  also  converted  from  microns  to  centimeters  to 
conform  with  other  test  data. 


4.2.5.  Subroutine  DATAIN 


This  subroutine  inputs  the  data  from  a  previously  generated  file 
into  a  dummy  array.  Under  the  control  of  subroutine  ORIENT,  the  data  are  then 
transferred  into  the  appropriate  3-dimensional  array. 

4.2.6.  Subroutine  DIFFER 

This  subroutine  provides  the  control  for  the  calculation  of  the 
first  and  second  derivatives.  The  order  of  processing  is: 

•  Displacement  is  inputted, 

•  Velocity  is  generated  using  an  FIR  differentiator, 

•  Velocity  is  smoothed  using  an  FIR  lowpass  filter, 

•  Acceleration  is  generated  using  an  FIR  differentiator,  and 


o 


Acceleration  is  smoothed  using  an  FIR  lowpass  filter. 


4.2.7. 


Subroutine  DISS4 


This  subroutine  generates  the  plot,  using  DISSPLA  (Version  8.2), 
and  stores  it  on  a  file  named  PLFILE.  This  plot  file  can  then  be  viewed  and 
printed  at  the  computer  terminal  and/or  sent  to  the  CALCOMP  plotter. 


4.2.8. 


Subroutine  FILT 


This  subroutine  inputs  the  FIR  coefficients  and  convolutes  the 
data  under  the  control  of  the  variable  I TYPE.  If  ITYPE=1,  the  data  are 
lowpass  or  highpass  filtered;  if  ITYPE=2,  the  data  are  differentiated.  Note 
that  the  processed  data  are  returned  in  the  same  array  in  which  they  were 
inputted. 


4.2.9. 


Subroutine  ORIENT 


This  subroutine  inputs  the  data  related  control  parameters, 
reads  NCH  sets  of  data  into  the  appropriate  part  of  the  3-dimensional  array 
DAT,  converts  the  data  into  engineering  units  if  ICVT  4  0,  and  lowpass  filters 
the  data  if  IFILT  4  0.  If  IDERIV  =  1,  displacement  data  are  inputted  and  the 
velocities  and  accelerations  are  calculated;  if  IDERIV  =  0,  all  three  are 
inputted.  For  each  data  set,  the  parameters  MIN,  NIN,  and  IU  are  read  in.  MIN 
and  NIN  are  indices  which  control  the  placement  of  the  data  set  into  the 
three-dimensional  array,  as  in  DAT ( MIN , NIN , L) ,  and  IU  indicates  on  which  tape 
unit  the  data  file  resides. 

4.2.10.  Subroutine  PL0T4 

This  subroutine  provides  the  control  over  what  is  plotted  by 
subroutine  DISS4.  Also,  the  data  must  he  put  into  1— dimensional  arrays  to  be 
plotted. 


4.2.11, 


Subroutine  WRDAT 


This  subroutine  writes  all  the  data  that  have  been  processed 
onto  a  file  on  tape  unit  18.  This  file  is  then  saved  and  used  as  an  input  to 
the  second  program. 

4.3.  Computer  Program  MUZPRED 

The  flowchart  for  MUZPRED  is  shown  in  Figure  4.2.  The  main  program 
inputs  the  control  variables  and  then  calls  the  subroutines.  The  control 
variables  are:  NCH  which  is  the  number  of  channels  to  be  inputted  (usually 
6),  NP1DTP,  NPLOTC  and  NPIDTU  which  contain  the  number  of  plots  to  be 
generated  in  subroutines  PREDCT,  CALCZ,  and  UNITV,  respectively,  and  IPR.INT 
which  controls  the  writing  of  an  output  data  file.  A  few  other  data  related 
variables  are  inputted:  LKEY ,  in  particular,  which  is  defined  as  the  index  in 
the  input  data  array  that  corresponds  to  the  time  when  the  midpoint  of  the 
rotating  band  of  the  projectile  passes  sensor  1.  When  plots  are  generated, 
time  is  zeroed  at  the  time  of  LKEY. 


17 


Figure  4.2.  Flowchart  For  Program  MUZPRED 


4.3.1.  Subroutine  GALCZ 

The  data  base  [5] [6]  on  which  the  development  of  this  analysis 
program  relied  did  not  contain  muzzle  recoil  motion.  Hence,  the  axial  muzzle 
motion  analysis  based  on  the  procedure  outlined  in  Eqs.  (3.10)  to  (3.12)  was 
not  incorporated  into  the  compute £  programs.  However,  since  this  motion  was 
obtained  by  an  earlier  experiment  using  the  same  gun  fixture,  projectile 
type  and  propellant  charge,  a  provision  has  been  added  to  this  subroutine  to 
scale  an  externally  supplied  function  for  the  muzzle  recoil  velocity  to  the 
projectile  in-bore  travel  and  gas  pressure  history.  This  function,  a  5-th 
order  polynomial,  was  generated  by  least  squares  model  fitting  of  the  average 
of  three  data  sets.  Its  coefficients  are  read  into  array  C.  The  times  of 
beginning  of  recoil  motion,  LZM,  maximum  pressure,  LPM,  and  shot  exit,  LEM,  of 
this  model  function  are  inputted  as  well  as  the  corresponding  last  two  of  the 
data  set,  LPD  and  LED.  After  proportionally  scaling  LZD ,  the  scaled  time 


Axial  muzzle  displacement  data  recorded  about  four  calibers 
back  of  muzzle  face  with  a  two-dimensional  displacement 
transducer  (ZIMMER)  were  graceously  supplied  to  BRL  by  Mr. 
T.O.  Andrews,  Royal  Armament  Research  and  Development 
Establishment,  Fort  Halstead,  Kent,  UK,  March  1985. 
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series  for  the  muzzle  recoil  velocity  is  computed  by  assigning  the  value  zero 
to  all  time  points  when  £<LZD  and  the  value  of  LED  to  all  time  points 
when  £>LED  .  The  displacement  array  is  then  obtained  by  numerical  integration 
of  the  velocity  over  time. 

The  distances  of  the  sensors  from  the  muzzle  are  read  into  a  3x2  array  Z 
along  with  the  reference  distance  ZREF.  At  each  point  in  time  a  least  squares 
fourth-degree  polynomial ,  with  the  second-  and  third-degree  coefficients  set 
equal  to  zero,  was  fit  to  the  data  of  all  three  sensors  according  to  Eqs. 

(3.4)  through  (3.9).  The  horizontal  and  vertical  fits  are  done  separately  and 
the  results  of  each  fit  are  evaluated  at  ZREF  and  saved  in  array  DAT. 

At  this  point  the  information  in  array  DAT  is: 


DAT( 1  ,1  ,  L) 

C0EF(3)  -  horizontal, 

DAT( 1  ,2  ,L) 

C0EF(3)  -  vertical. 

DAT( 2 , 1  ,  L) 

displacement  at 

ZREF  -  horizontal, 

DAT(2,2,L) 

displacement  at 

ZREF  -  vertical, 

DAT( 3,1  ,L) 

slope  at  ZREF  - 

horizontal , 

DAT(3  ,2  ,L) 

slope  at  ZREF  - 

vertical , 

DAT( 4 ,1  ,L) 

C0EF(1)  -  horizontal, 

DAT(4,2,L) 

C0EF(1)  -  vertical. 

DAT(5  ,1,L) 

C0EF(2)  ~  horizontal, 

DAT(5,2,L) 

C0EF(2)  -  vertical, 

DAT(6 ,  1  ,L) 

recoil  velocity. 

,  and 

DAT(6,2,L) 

recoil  displacement . 

Next,  subroutine 

UNITV  is  called  to 

calculate  the  unit 

discrete  time,  using  the  slopes.  That  procedure  is  described  in  Section 
4.3.11.  During  this  procedure,  the  arrays  DAT(6,1,L)  and  DAT(6,2,L)  are 
overwritten. 


The  displacement  at  ZREF  is  then  differentiated  to  form  the  velocity, 
lowpass  filtered  and  saved  also  in  array  DAT,  overwriting  DAT(3,1,L)  and 
DAT(3,2,L).  At  this  point  the  information  in  array  DAT  is: 


DAT(1 ,1 ,L) 
DAT( 1 ,2 ,L) 
DAT( 2 , 1 ,L) 
DAT( 2  ,2  ,L) 


C0EF(3)  -  horizontal, 

C0EF(3)  -  vertical, 
displacement  at  ZREF  -  horizontal, 
displacement  at  ZREF  -  vertical , 


19 


DAT(3 ,1 ,L) 

velocity  at  ZREF  -  horizontal. 

DAT( 3 ,2 ,L) 

velocity  at  ZREF  -  vertical, 

DAT( 4 ,1  ,  L) 

C0EF(1)  -  horizontal, 

DAT( 4  ,2  ,  L) 

C0EF(1)  -  vertical, 

DAT(5 ,1 ,L) 

C0EF(2)  -  horizontal, 

DAT(5 ,2 ,L) 

C0EF(2)  -  vertical, 

DAT (6,1 , L) 

e31  ’  and 

DAT(6,2,L) 

e32  * 

The  plots 
and  NDAT .  Any 

are  then  generated  as  specified  by  NPLOTC  and 
inputs  required  for  plotting  reside  on  a  file 

the  arrays  MDAT 
on  tape  unit  7 • 

4.3.2. 

Subroutine  C0EF1 

parameters  of 

This  subroutine  inputs  the  coefficient  array  and  other  necessary 
a  previously  generated  FIR  lowpass  or  highpass  filter. 

4.3.3. 

Subroutine  C0EF2 

parameters  of 

This  subroutine  inputs  the  coefficient  array  and  other  necessary 
a  previously  generated  FIR  differentiator. 

4.3.4. 

Subroutine  CONVOL 

This  subroutine  does  a  linear  convolution  of  the  NFILT 
coefficients  with  a  data  array  of  length  NFILT.  This  process  then  proceeds  in 
moving  arc  fashion  along  the  entire  data  array.  The  convolution  at  the 
endpoints  are  accounted  for  by  reflection  of  the  data  at  the  endpoints.  See 
the  listing  in  Appendix  C  for  more  details. 

4.3.5. 

Subroutine  DATAIN 

This  subroutine  inputs  the  data  from  program  MUZM040  into  a 
three-dimensional  array  as  specified  by  variables  M  and  N.  The  input  data 
array  dimension  is  4100,  but,  because  of  central  memory  size  limitations,  the 
data-processing  array  DAT  is  restricted  to  2500  points. 

4.3.6.  Subroutine  DATA4 

This  subroutine  prepares  the  data  arrays  for  plotting.  Special 
provision  is  made  to  zero  the  data  at  the  time  LKEY  mentioned  in  the  main 
program;  the  abscissa  is  zeroed  on  the  time  plots  and  both  abscissa  and 
ordinate  on  the  x-y  plots. 
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4.3.7 


Subroutine  DIFFER 


This  subroutine  calls  the  subroutines  required  for 
differentiating  if  IUD  4  0  and  then  calls  the  ones  required  for  lowpass  or 
highpass  filtering  if  IUF  4  0. 

4.3.8.  Subroutine  DISS4 

This  subroutine  generates  the  plot,  using  DISSPLA  (Version  9.0), 
and  stores  it  on  a  local  file  named  META.  Up  to  three  curves  can  be  plotted 
and  options  are  available  for  color,  scales  and  labels.  For  more  details  see 
the  subroutine  listing. 

4.3.9.  Subroutine  FILT 

This  subroutine  reads  in  the  FIR  lowpass  or  highpass  filter 
coefficients  and  calls  the  subroutine  CONVOL  to  convolute  the  data.  Note  that 
the  processed  data  are  returned  in  the  same  array  in  which  they  were  inputted. 

4.3.10.  Subroutine  PREDCT 


This  subroutine  does  a  straight  line  extrapolation  from  the  data 
of  sensors  2  and  3  to  the  position  of  sensor  1,  the  variables  DZ31  and  DZ32 
being  the  distances  from  sensor  3  to  sensor  1  and  2,  respectively.  The  values 
of  the  prediction  can  then  be  compared  with  the  ones  of  sensor  1;  the 
magnitudes  and  phase  angles  of  both  data  sets  are  calculated  to  aid  in  this 
comparison. 

In  the  vicinity  of  LKEY  the  displacements  are  averaged  over  short  time 
spans  and  their  magnitudes  and  phase  angles  are  calculated,  stored,  and 
printed. 

At  this  point  the  information  in  array  DAT  is: 


DAT( 1 ,1 ,L) 
DAT( 1 ,2 ,L) 
DAT( 2 ,1 ,L) 
DAT( 2 ,2 ,L) 
DAT( 3 , 1 ,L) 
DAT( 3  ,2  ,L) 
DAT( 4 , 1  ,L) 
DAT( 4  ,2  ,L) 
DAT( 5 , 1  ,L) 
DAT(5  ,2  ,L) 


sensor  1  displacement  -  horizontal, 
sensor  1  displacement  -  vertical , 
sensor  2  displacement  -  horizontal, 
sensor  2  displacement  -  vertical  , 
sensor  3  displacement  -  horizontal, 
sensor  3  displacement  -  vertical, 
prediction  displacement  -  horizontal, 
prediction  displacement  -  vertical, 
prediction  magnitude, 
prediction  phase  angle,  degrees, 
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DAT(6 ,1 ,L) 
DAT(6 ,2 ,L) 


sensor  1  magnitude,  and 


sensor  1  phase  angle,  degrees. 

NPLOTP  plots  are  then  generated  with  the  input  arrays  MDAT  and  NDAT 
controlling  what  is  plotted.  Any  inputs  required  for  plotting  reside  on  a 
file  on  tape  unit  3. 

4.3.11.  Subroutine  UNITV 

This  subroutine  calculates  the  translation  and  rotation 
parameters  as  described  in  Section  3.2.3,  Eq.  (3.16)  in  particular  .  The 
three-dimensional  array  DAT,  at  entry  to  the  subroutine,  contains  the 
information  as  shown  in  subroutine  CALCZ.  The  array  is  overwritten  again  in 
this  subroutine  as  follows: 


DAT(1 ,1  ,L) 

e!3  ,  j  =  1,2, 

or  3 , 

DAT( 1  ,2  ,L) 

( irrelevant) , 

DAT( 2 , 1 ,  L) 

(irrelevant) , 

DAT(2 ,2 ,L) 

( irrelevant) , 

DAT(3 ,1  ,L) 

slope  at  ZREF-horizontal , 

DAT(3 ,2 ,L) 

slope  at  ZREF-vertical , 

DAT (4,1 , L) 

(irrelevant) , 

DAT(4,2,L) 

( irrelevant) , 

DAT(5,1.L) 

(irrelevant) , 

DAT(5,2,L) 

(irrelevant) , 

DAT(6 ,1 ,L) 

eji  >  j  -  1.2, 

or  3 ,  and 

DAT(6 ,2 ,L) 

ej2  ,  j  =  1,2, 

or  3* 

NPLOTU  plots  are  then  generated  as  specified  by  the  input  arrays  MDAT  and 
NDAT.  Any  inputs  required  for  plotting  reside  on  a  file  on  tape  unit  10. 

4.3.12.  Subroutine  WRFILE 

This  subroutine  takes  data  from  a  three-dimensional  array  and 
transfers  all  or  part  of  one  column  into  a  one— dimensional  array  which  is  then 
written  on  an  output  file  located  on  tape  unit  9.  This  procedure  is  done  NARR 
times  and  the  columns  to  be  written  are  chosen  according  to  the  input  arrays 
MDAT  and  NDAT.  LST  and  LSP  are  the  starting  and  stopping  indices, 
respectively,  of  the  data  to  be  written. 
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5.  CONCLUSIONS 


The  methodology  described  in  the  report  allows  the  exploitation  of  the 
SDT  technology  for  the  determination  of  the  muzzle  motion  during  the  in-bore 
and  launch  motion  of  the  projectile  with  a  very  high  accuracy.  The 
theoretical  foundation  is  general.  It  contains  all  the  equations  and 
procedures  necessary  for  the  analysis  of  three-dimensional  time  series 
presenting  the  tube  displacement  in  the  longitudinal  as  well  as  in  the 
horizontal  and  vertical  directions  from  SDT  type  instrumentation  at  k- 
locations.  This  allows  the  determination  of  five  of  the  six  degrees-of- 
freedom  parameters  as  a  function  of  time  for  all  points  lying  on  the  tube  axis 
in  the  spatial  measurement  domain.  Only  the  rotation  about  the  longitudinal 
tube  axis  is  not  obtainable.  For  most  gun  accuracy  related  investigations 
however,  the  knowledge  of  this  parameter  is  not  required,  since  one  is  mainly 
interested  in  the  transverse  velocity  and  the  aiming  direction  of  the  muzzle 
and  in  their  temporal  derivatives  during  shot  exit. 

Currently  the  computer  program  is  explicitly  set  up  for  the  analysis  of 
transverse  data  collected  at  three  locations  at  the  muzzle  end  of  the  tube 
while  allowing  for  the  muzzle  recoil.  It  can  easily  be  extended,  if 
necessary,  to  n  >  3  locations  and  to  include  determination  of  the  axial  muzzle 
motion  from  respective  SDT  data,  employing  the  outlined  theoretical 
formulation . 
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APPENDIX  A 


BRIEF  REVIEW  OF  CANNON  MOTION  MEASUREMENT  TECHNIQUES 
A.  1  Introduction 

Various  instrumentation  and  measurement  techniques  are  available  to 
monitor  the  gun  tube  motion  during  the  in-bore  travel  and  launch  of  a 
projectile.  Most  methods  are  based  on  the  utilization  of  streak  photography, 
electroopt ical  displacement  transducers,  optical  levers,  electromagnetic 
proximity  transducers,  interferometry,  electromechanical  devices,  strain 
gauges,  accelerometers,  etc.  The  authors,  not  being  aware  of  any  publication 
containing  a  review  of  experimental  methods  for  measuring  tube  motion,  felt 
compelled  for  reasons  of  completeness  to  provide  a  very  brief  description  of 
the  above  mentioned  methods  with  some  references  to  applications, 

A, 2.  Streak  Photography 

Streak  photography  [A-l ]  is  probably  one  of  the  best  photographic 
techniques  for  recording  the  motion  history  of  a  vibrating  tube.  In  this 
instrumentation,  the  image  formed  by  a  small  reflector  bonded  to  the  surface 
is  recorded  on  film  by  a  moving-film,  shutterless  camera  in  the  form  of  a 
continuous,  wavy  line  against  a  contrasting  background.  This  time- 
displacement  record  can  be  evaluated  with  an  accuracy  of  the  order  of  10  cm 
and  a  time  resolution  of  about  1  msec.  Due  to  the  extra  labor  and  time 
involved  in  processing  the  optical  data  this  technique  is  relatively  costly 
over  techniques  which  utilize  electronic  recording. 

A. 3*  Electrooptical  Displacement  Transducers 

Electroopt ical  displacement  transducer  instrumentation  combines  the 
advantage  of  optics  and  electronic  recording  by  incorporating  an  optical 
system  in  an  electronic  circuit.  Two  different  arrangements  are  described 
here. 


The  first  one  employs  a  light-dark  discontinuity  on  the  cannon  surface 
which  is  imaged  by  a  lens  system  onto  a  photocathode  which  converts  the 
optical  image  into  an  electron  beam.  This  beam,  in  turn,  is  imaged  via  an 
electromagnetic  field  onto  the  center  of  a  small  aperture  and  thereafter 
amplified  by  an  electron  multiplier.  The  resulting  current  is  used  as  a  servo 
loop  control  for  the  electromagnetic  field  generator.  The  correcting  current, 
necessary  to  keep  the  electron  image  centered  on  the  aperture,  is  recorded  as 
a  measure  of  target  deflection.  Displacement  range  and  resolution  are  a 
function  of  the  focal  length  of  the  lens  system  and  its  distance  from  the 
target.  In  the  US,  one-dimensional  electrooptical  transducers  manufactured  by 
Optron,  a  Division  of  Universal  Technology  Inc.,  of  Woodbridge,  Connecticut 
have  been  successfully  employed  for  measuring  gun  dynamics  [A-2 ] [A-3 ] [A-4 3 • 

In  Western  Europe  two-dimensional  electrooptical  transducers  manufactured  by 
Zimmer  OHG,  Darmstadt,  FRG,  are  employed  for  the  same  purpose  [A-5 ] [A-6 ] [A-7 ] 
and  as  a  standard  for  calibrating  other  less  expensive  instrumentation  for 
measuring  tube  motion. 

In  the  second  arrangement,  a  collimated  light  beam  passes  over  a  knife 
edge  and  partially  illuminates  a  large  area  light-sensitive  photodiode  through 
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a  narrow  slit.  The  knife  edge,  e.g.,  a  razor  blade,  is  mounted  to  the  gun 
tube,  parallel  to  the  tube  axis  and  perpendicular  to  the  light  beam,  and  is 
positioned  such  that  it  partially  blocks  the  light  beam  passing  through  the 

slit.  For  instance,  if  the  knife  edge  is  arranged  horizontally  in  the 

vertical  plane  of  the  tube  axis,  any  vertical  tube  motion  will  move  the  edge 
up  or  down  and  regulate  the  light  striking  the  sensor.  The  detected  signal  is 

then  amplified  and  recorded.  The  signal  change  is  directly  proportional  to 

the  gun  tube  displacement  in  the  vertical  direction.  By  using  two  light 
sources  and  corresponding  sensors  positioned  orthogonally,  both  the  horizontal 
and  vertical  components  of  the  tube  motion  can  be  measured  [A— 8 ] [A— 9 ] • 

A.4.  Optical  Lever  Instrumentation 

Optical  lever  instrumentation  is  ideally  suited  for  monitoring  local 
angular  deflection  of  gun  tubes  [A-10].  It  has  recently  been  integrated  with 
electronic  recording  [A— 1 1 ]  and  applied  to  the  measurement  of  the  muzzle- 
pointing  direction  [A— 12] [A— 13] [A— 14].  In  this  measurement  technique  a 
collimated  light  beam  is  incident  on  and  reflected  by  a  mirror  rigidly 
attached  to  the  gun  tube  and  imaged  via  a  lens  system  onto  a  two-dimensional, 
position-sensing  photodetector.  As  the  gun  tube  changes  its  curvature,  the 
mirror  rotates  with  it,  thereby  changing  the  direction  of  the  reflected  light 
beam  and,  thus,  the  location  of  the  light  spot  on  the  photodetector.  The  x 
and  y  components  of  the  induced  photocurrent  are  recorded  as  a  function  of 
time  and  converted  into  a  displacement  history  with  respect  to  the  electric 
center  of  the  sensor.  Superposition  of  the  geometry  of  the  optical  lever 
setup  yields  the  local  angular  deflection  history.  This  type  of 
instrumentation  is  employed  in  muzzle  reference  systems  for  tank  guns. 

A. 5m  Electromagnetic  Proximity  Transducers 

Electromagnetic  proximity  transducers  are  increasingly  being  used  for 
measuring  tube  displacements  [ A— 1 5 ] [A-l 6 ] [A-l 7 ] .  Most  operate  on  the  eddy 
current  principle:  a  varying  electromagnetic  field,  usually  in  the  radio 
frequency  range,  is  generated  by  an  oscillator  and  radiated  from  its  active 
inductance  coil/antenna  to  a  proximate  surface  location  or  protrusion  of  the 
tube.  Eddy  currents  are  induced  in  the  metallic  surface  and  generate  an 
electromagnetic  field  which  then  couples  back  into  the  active  inductance 
coil.  As  the  surface  vibrates,  the  recoupling  to  the  inductance  coil  changes, 
which,  in  turn,  causes  a  change  in  the  oscillator  impedance  and,  thus,  a 
modulation  in  the  oscillator  current  frequency  and/or  current  amplitude.  This 
modulation  is  recorded  as  a  function  of  time  and  transferred  into  a 
displacement  history.  Inductive  proximity  probes  which  can  readily  be  adapted 
to  the  measurement  of  gun  dynamics  are  commercially  available. 

The  Schmidt  displacement  transducer  (SDT)  described  in  this  report  is  a 
proximity  transducer,  which  integrates  over  the  tube  perimeter,  thus 
suppressing  contributions  from  local  vibrations  and  yielding  a  large  signal- 
to-noise  ratio. 

A.  6.  Interferometric  Measurement  Techniques 

Interferometric  measurement  using  laser  as  well  as  microwave  techniques 
can  also  be  used  to  record  certain  dynamic  parameters  of  tube  motion.  In  this 
technique  a  monochromatic  coherent  electromagnetic  wave  is  divided  by  a  beam 
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splitter  into  two  parts  which  travel  different  paths  and  recombine  to  form 
interference  fringes.  If  one  of  the  beams  is  reflected  by  a  moving  object, 
which  in  our  case  is  the  tube  surface,  well-defined  maxima  and  minima  in  the 
fringe  intensity  are  produced  and  recorded  either  photographically  or 
electrically.  MW  interferometry  represents  a  convenient  tool  to  monitor  the 
axial  recoil  motion  history  of  guns  at  the  breech  as  well  as  at  the  muzzle. 
Using  an  appropriate  experimental  setup,  the  latter  event  can  even  be  obtained 
together  with  projectile  travel  from  in-bore  MW  interferometry  [A-18].  In  the 
visible  spectrum,  D.  Warken  [ A— 1 9 ] [A-20]  pioneered  the  application  of 
holographic  interferometry  and  laser  speckle  photography  for  the  investigation 
of  the  spatial  displacement  field  of  local  tube  surfaces  in  clearly  defined 
time  intervals.  Though  the  experimental  arrangement  for  providing  a  time 
sequence  of  holograms  or  speckle  photographs  may  be  quite  cumbersome  to  set  up 
and  the  recorded  information  difficult  to  analyze,  these  optical  methods  are 
well  suited  for  the  detailed  investigation  of  the  displacement  field  of  those 
locations  on  the  gun  surface  which  have  inhomogeneities  in  geometry  or  in 
material  properties.  For  typical  gun  tube  vibration  measurements  however, 
they  are  too  complex  and  too  detailed.  On  the  contrary,  laser  interferometry 
with  electronic  recording  is  very  attractive  for  the  temporal  measurement  of 
local  tube  motion.  A  multitude  of  interferometer  configurations  have  been 
developed  to  meet  numerous  experimental  requirements.  Three  of  them, 
generally  identified  by  the  acronyms  DISAR  [A-21],  VISAR  [A-22]  and  TRANSAR 
[A-23],  should  easily  be  adaptable  for  the  measurement  of  gun  motion  [A-24]. 
The  first  is  basically  a  Michelson  interferometer  in  which  the  light  source  is 
a  high-powered  laser  and  the  target  a  diffuse  scatterer.  The  backscattered 
light  is  collected  by  a  telescope  and  brought  to  interference  with  the 
unscattered  reference  beam.  In  the  second  configuration,  essentially  a  laser 
Doppler  velocimeter,  the  backscattered  light  is  again  collected  and  collimated 
by  a  lens  system  and  is  divided  into  two  beams.  After  sending  one  beam  around 
an  optical  delay  leg  it  is  recombined  with  the  undelayed  beam.  The 
interference  fringe  history  thus  formed  is  proportional  to  the  velocity  change 
of  the  reflecting  surface.  In  the  third  configuration,  the  TRANSAR,  the 
backscattered  light  is  collected  at  two  different  directions,  usually  in  a 
diagonally  opposite  setup,  via  telescopes  and  optically  heterodyned  to  yield 
the  temporal  displacement/velocity  of  the  scattering  surface  location 
orthogonal  to  the  direction  of  the  incident  beam  and  in  the  plane  formed  by 
the  direction  in  which  the  scattering  is  observed. 

A.  7.  Electromechanical  Devices 

A  recent  innovation  of  an  electromechanical  device  is  the  "tuning  fork" 
transducer  developed  by  S&D  Dynamics  [A-25].  In  this  measurement  technique,  a 
cylindrical  section  of  the  cannon  is  sandwiched  between  two  prongs  of  a  tuning 
fork-like  mechanical  device  with  its  stem  rigidly  attached  to  a  nonmovable 
mass.  As  the  tube  moves  towards  one  of  the  prongs,  it  forces  the  prongs  to 
move  with  it,  which,  in  turn,  causes  a  bending  of  the  stem.  This  deformation 
is  then  picked  up  by  strain  gauges  mounted  on  the  stem  and  recorded. 

A. 8.  Strain  Gauge  Instrumentation 

Strain  gauge  instrumentation  usually  employs  a  matrix  of  temperature- 
compensating  strain  transducers  appropriately  aligned  and  mounted  on  the  outer 
surface  of  the  cannon  to  measure  local  tube  strains  [A-6][A-26].  Each  strain 
gauge  forms  an  active  leg  in  a  Wheatstone  bridge.  The  strain  contribution  to 
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tube  bending,  however,  must  be  separated  analytically  from  that  produced  by 
the  gas  pressure,  the  gun  recoil,  and  the  dynamic  contact  of  the  projectile 
with  the  bore  surface  during  projectile  in-bore  travel  and  launch,  before  it 
can  be  translated  into  tube  curvature. 

A. 9.  Accelerometer  Instrumentation 

The  accelerometer  instrumentation  as  customarily  employed  [A-17 ] [ A-26 ] [ A- 
27][A-28]  uses  piezoelectric  accelerometers  which  are  bonded  to  the  outer 
surface  of  the  cannon  in  a  matrix  arrangement  to  record  local  tube 
accelerations.  Because  of  the  unfortunate  susceptibility  of  piezoelectric 
transducers  to  high  frequency  shocks  which  can  produce  temporary  step  function 
like  zero  shifts  in  the  signals,  their  application  for  monitoring  cannon 
vibration  is  limited. 
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ESTABLISHMENT  OF  GANNON  FLEXURE  WITH  RESPECT  TO  BORE  SIGHT  REFERENCE  LINE 
PRIOR  TO  SHOT  START 

There  are  a  few  gun-related  observables  which  should  or  must  be 
determined  just  prior  to  the  commencement  of  fire.  They  represent  the  initial 
values  of  gun  flexure  and  aiming  point.  If  their  determination  is 
simultaneously  done  with  the  measurement  of  the  tube  displacements,  we  can 
establish  a  frame  of  reference  where  these  three  sets  of  observables  are 
correlated.  Hence,  it  is  presumed  that  a  muzzle/ tube  displacement  measurement 
arrangement  as  described  in  Section  3.1  is  an  integral  part  of  the  overall 
measurement  setup  and  is  used  to  monitor  the  tube  flexure  during  the 
establishment  of  the  aiming  point  and  the  initial  tube  curvature. 

B.l.  Gun  Aiming  Point/ System  Frame  Of  Reference 

In  many  gun  accuracy-related  investigations,  we  are  interested  in  the 
relationship  between  the  projectile  impact  and  the  initial  aiming  point  of  the 
gun.  Restricting  the  discussion  to  direct  fire  weapons,  we  can  employ  a 
boresight  or  an  alignment  telescope  for  the  establishment  of  the  aiming  point 
or  points  on  the  target.  Because  of  the  simplicity  of  the  measurement,  the 
intersection  of  the  line  of  sight  given  by  a  boresight  placed  into  the  muzzle 
or  breech  section  of  the  cannon  with  the  target  witness  board  is  habitually 
used  as  the  reference  point.  However,  any  other  section  of  the  bore  may  be 
used  for  establishing  the  line  of  aim.  For  example,  the  line  of  sight  which 
passes  through  the  center  of  the  breech  and  the  center  of  the  muzzle  is  also 
employed  as  a  line  of  aim  as  reference. 

In  the  following  we  assume  that  the  line  of  aim  is  established  by  a 
boresight  placed  into  the  muzzle  and  presents  the  mean  of  many  measurements  in 
which  the  boresight  has  systematically  been  rotated  and,  if  necessary, 
reseated.  By  statistically  averaging  the  observations,  we  can  practically 
eliminate  gun  locality  and  human  related  bias.  At  the  same  time  the  aiming 
points  are  recorded,  we  may  also  record  the  muzzle  displacements  and  calculate 
from  their  averages  the  location  and  the  slope  of  the  muzzle.  These 
expressions  are  respectively, 


SA(v)  -  S(v  ,z=0,A  ),  v=l  ,2  and  e^  =  e»(X),  i=l,2,3 


(B.l) 


using  the  algorithm  described  in  Section  3.2.  Because  of  the  free  beam 
condition  at  the  muzzle  face,  the  muzzle  end  is  basically  a  rigid  hollow 
cylinder.  This  allows  us  to  postulate  a  frame  of  reference  which  has  its 
origin  at  (SA(1),  SA(2),0)  and  its  third  axis  pointed  to  the  aim  point.  It 
has  to  be  noted  that  this  line  of  aim  is  not  the  actual  one  at  shot  start, 
since  it  has  been  established  with  the  additional  weight  of  the  boresight  in 
the  muzzle.  If  the  boresight  is  removed,  the  tube  flexure  will  readjust 
itself  to  a  new  slightly  different  stress  equilibrium,  thereby  changing  its 
muzzle  shape.  However,  with  the  system  frame  of  reference  defined,  any  change 
in  the  muzzle  motion  can  be  expressed  relative  to  it. 
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B.2.  Bore  Straightness /Axis 

Generally,  gun  tubes  are  not  straight*  Reasons  for  this  include  residual 
stresses,  variation  in  wall  thickness,  machining  accuracy  and  tolerances, 
eccentricity  of  rifling,  storage  position,  environment,  etc.  Also  large 
thermal  distortion  can  be  induced  by  solar  irradiation.  The  sun-exposed  side 
of  the  cannon  will  warm  up  and  expand  with  respect  to  the  shadowed  part.  This 
produces,  in  addition  to  the  natural  crookedness  and  the  gravitational  droop, 
a  curvature  in  the  cannon  with  its  center  of  curvature  lying  in  the  shadowed 
plane,  thus  bending  the  muzzle  away  from  the  sun.  Hence,  it  is  important  to 
measure  the  initial  conditions  of  the  tube  bending  immediately  prior  to  the 
commencement  of  fire  for  gun  accuracy,  diagnostic,  and  model  simulation 
purposes • 

Up  till  now  tube  straightness  or  bend  has  been  measured  with  an  alignment 
telescope  and  an  illuminated  bull's-eye  target  which  is  moved  through  the  tube 
[B-l].  The  telescope  is  placed  on  an  adjustable  mount  near  the  muzzle  or 
breech  of  the  tube  and  adjusted  so  that  the  line  of  sight  passes  through  the 
center  of  the  bore  at  the  commencement  of  the  rifling  or  forcing  cone  and  the 
muzzle.  The  deviations  from  the  line  of  sight  are  measured  at  various 
positions  along  the  tube  axis  in  both  horizontal  and  vertical  directions  with 
an  accuracy  of  about  .025mm  by  means  of  micrometer  attachments  on  the 
telescope.  This  cumbersome  optical  measurement  technique  is  amenable  to 
automatization  by  introducing  a  laser  beam  as  line  of  reference,  replacing  the 
bull's  eye  target  by  a  two-dimensional  electrooptical  position  sensor, 
mechanizing  the  push  or  pull  of  the  target  carriage  through  the  bore,  and 
digitally  recording  the  beam  location  as  a  function  of  the  axial  displacement 
of  the  target.  Such  a  device  has  recently  been  introduced  by  Watervliet 
Arsenal  to  measure  the  tube  curvature  during  and  after  manufacturing.  The 
target  carriage  is  centered  in  the  bore  by  mechanical  springs  which  limits  the 
accuracy  of  the  measurement.  Because  gun  dynamics  model  simulations  strongly 
indicate  that  muzzle  motion  at  shot  exit  is  very  sensitive  to  the  tube  bend  at 
shot  start,  the  inaccuracy  introduced  by  the  mechanical  alignment  mechanism  of 
the  target  carriage  with  the  bore  center  must  be  reduced  as  much  as  possible. 

For  smoothbore  guns,  this  can  be  done  with  relative  ease  by  mounting  a 
SDT  system  into  the  outer  surface  of  the  probe  concentric  and  coplanar  with 
the  photodiode  and  recording  the  respective  signals  concurrently  as  a  function 
of  the  axial  displacement.  The  SDT  instrumentation  provides  the  location  of 
the  geometric  center  of  the  bore  relative  to  the  probe  center  and  the  position 
sensor  gives  the  location  of  the  reference  laser  beams  with  respect  to  the 
probe  center.  By  introducing  a  parallel  plate  which  is  mounted  on  a  motor- 
driven,  computer-controlled  goniometric  cradle  or  rotational  displacement 
device  into  the  light  path,  we  can  displace  the  light  beam  in  the  plane  of 
rotation  in  a  controlled  way.  By  letting  the  parallel  plate  oscillate  about 
its  zero  position  normal  to  the  beam  and  recording  the  beam  location  at  the 
zero  transition  and  the  extreme  displacements,  we  even  can  account  for  any 
rotational  movement  of  the  probe  about  the  axis  and  uniquely  correlate  the 
bore  axis  relative  to  the  reference  coordinate  system  given  by  the  laser  beam 
arrangement  (Figure  B.l).  This  measurement  will  yield  the  geometric  center  of 
the  bore,  a,  with  regard  to  the  reference  beam  and  as  a  function  of  the  axial 
position  Q: 
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a(0  =  (5,n,5) 


(B.  2) 


CENTER 
OF  PROBE 


Figure  B.l.  Accurate  Measurement  Of  Bore  Curvature  Prior  To  Shot  Start 


Measuring  the  tube  flexure  at  the  same  time  yields 


S*(v,z,£)  ,  v=l ,2, 


(B.3) 


where  z  and  4  are  the  locations  of  the  SDT  instrumentation  and  the  probe, 
respectively*  We  can  calibrate  the  tube  axis  as  defined  by  the  individual  SDT 
measurements  to  the  bore  axis  as 


S*(v,z,£)  =  S*(v,z3s)  -  [S*(v,z,z)  -  a^( z)],  v=l,2* 


(B*4) 


If  we  align  the  laser  reference  beam  collinear  with  the  line  of  aim  we  can 
avoid  minute  translation  and  rotation  in  correlating  the  two  reference  systems 
and  can  define  the  alignment  function  So(v,y)  appearing  in  Eq*  (3*2)  as 


So(v5y)  =  S*(v,Zy,zy)  -  av(zy) 


v— 1 j  2  j  P“1 $  2  y • • .m  • 


(B.5) 
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APPENDIX  C 
COMPUTER  PROGRAMS 
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APPENDIX  C 


COMPUTER  PROGRAMS 


C.l.  Program  MDZM040 

C.1.1.  Listing  Of  Job  Control  Language  Of  MFA  File  MUZMO40/UN=BOOTS 

BOOTS, STMFZ,P5,T50,MS300000. 

ACCOUNT, PDXXX. 

REQUEST, PLFILE,*PF. 

REQUEST, TAPE7,*PF. 

REQUEST, TAPE 1 8 , *PF . 

ATTACH, DISSPLA,ID=DISSPLA. 

ATTACH, COMPRS,ID=DISSPLA. 

LIBRARY, *,DISSPLA,COMPRS. 

FTN,LCM=I,L=0. 

BEGIN, GETMFATI,  FILE, LF=TAPE1 , PF=BR05 12A ,UN=BOOTS . 

BEGIN , GETMFAU , FILE , LF=TAPE 2 , PF=BR05 1 1 A , UN=BOOTS . 

BEGIN, GETMFAU, FILE, LF=TAPE8,PF=BR0509A,UN=BOOTS. 

BEGIN, GETMFAU, FILE , LF=TAPE9 , PF=BR0508A, UN=BOOTS . 

BEGIN, GETMFAU , FILE , LF=TAPE 1 0 , PF=BR05 1 OA , UN=BOOTS . 

BEGIN , GETMFAU , FILE , LF=TAPE 1 1 , PF=BR0507A , UN=BOOTS . 

BEGIN, GETMFAU, FILE ,LF=TAPE4 , PF=DIFF1 , UN=BOOTS . 

BEGIN, GETMFAU, FILE, LF=TAPE17,PF=LPF62,UN=BOOTS. 

MAP, OFF. 

LGO. 

BEGIN, SAVMFAU, FILE, LF=TAPE18,PF=DVA05U,UN=BOOTS. 


APPENDIX  C 


C.1.2.  Listing  Of  FORTRAN  Program 
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10 


15 
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PROGRAM  MUZZLE ( INPUT , OUTPUT , TAPE 5=INPUT , TAPF6=0UTPUT , TAPE  1 , TAPE 1 7 

*  , PLF I LE , TAPE  2 , TAPE  3 , TA PE4 , TAPE 7 , TAPE R , TAPE  9 , TA PE 1 0 , TAPE  1 1 , TA PE 1 R ) 
C 

C  ANALYSIS  OF  MUZZLE  MOTION  DATA  -  PROGRAM  1 

C 

COMMON/ I /MIN ( 6 ) , NIU( 6 ) , NTS ( 24 ) , NCH , IDERIV , IFILT , NPT , LDEL 

*  , NLSQ, LABEL C 10) ,LABL( 3, 2, 24) , LCNT( 2,7) 

REAL  LDEL 

COMMON  DUM(l), 

*  X(5000) ,Y1 (5000) ,Y2(5000) ,Y3(5000) ,Y4(5000) ,Y5(5000) ,Y6(5000) 
COMMON/ LE VI /DTIME ,CT( 3, 2), CT 1(3, 2) ,CT2( 3,2) ,CT3(3,2) ,CVT(3,6) 

*  , TSTART ,  ZZ( 3 , 2 ) 

COMMON/LEV2C/DAT(3 , 2 , 5000) ,DATD(3 , 2 , 5000) ,DATDD(3 , 2 , 5000) 

LEVEL  2 , DAT , DATD , DATDD 
LABEL ( 1 )=10H  BOOTS 
LABEL(2)=10HBLDG. 390 
LABEL(3)=10HX  6121 
LABEL(4)=10HMUZZLE 
CALL  COMPRS 
CALL  SETDEV(0 , 6 ) 

READ(5,2)IPLT1,IPLT2, IPLT3.IPRT 
WRITE ( 6,1) IPLT1 , IPLT2 , IPLT3 , IPRT 
CALL  ORIENT 

IF(IPLT1.E0.1)CALL  PLOT4(DAT) 

IF( IDERIV .EO.O)GO  TO  100 
CALL  DIFFER 

IF( IPLT2 .EO. 1 ) CALL  PLOT4(DATD) 

IF(IPLT3.E0.1)CAI.L  PLOT4(DATDD) 

100  IF( IPRT.EO. 1 )CALL  WRDAT 
CALL  DONE PL 
STOP 

1  F0RMAT(1H1,'  IPLT1  =',I3,’  IPLT2  =',I3,'  IPLT3  IPRT  =',I3) 

2  F0RMAT(10I3) 

END 
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SUBROUTINE  C0EF1 ( NFILT,H, IU) 

READS  IN  COEFFICIENTS  OF  AN  FIR  LOWPASS  OR  HIGHPASS  FILTER 

DIMENSION  H(512) 

READ( IU) NFILT , FP ,TBWID , DP ,DS 
IF(lSW.NE.l ) WRITE (6 , 1 ) 

IF(ISW.NE.l ) WRITE (6 , 2 )NFILT, FP , TBWID , DP ,DS 

N=( NFILT+1 ) /2 

READ ( IU ) ( H( I ) , 1= 1 , N) 

IF( ISW.NE. 1 ) WRITE (6 ,3) 

IF( ISW.EQ.l )GO  TO  200 
DO  100  1=1 ,N,  5 
IST=I+4 

IF( 1ST .GT.N) IST=N 

WRITE (6 ,4)I,(H(II), 11= I, 1ST) 

100  CONTINUE 

WRITE(6 ,5) 

200  ISW=1 
RETURN 

1  F0RMAT( // , '  FILTER  PARAMETERS, ') 

2  FORMAT(5X, '  NFILT  =  ',15,'  FP  =  ',F10.5,'  TBWID  =  ',F10.5, 
*  *  DP  =  ' ,F10 .5  , '  DS  =  1 ,F10 .5) 

3  FORMAT(5X, '  FILTER  COEFFICIENTS') 

4  FORMAT(5X,I5,5E15.8) 

5  FORMAT(//) 

END 
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SUBROUTINE  COEF 2 ( NF ILT , H , IU , JTY PE ) 


C 


READ  IN  COEFFICIENTS  OF  AN  FIR  DIFFERENTIATOR 


DIMENSION  H(A4) ,EDGE( 20) , FX( 10) ,WTX(10) 

IF( ISW.NE . 1 ) WRITE ( 6 , 1 ) 

READ (IU)NFILT, JTYPE , NBANDS 

IF( ISW. NE.1)WRITE(6, 2 )NFILT, JTYPE, NBANDS 

JB=2*NBANDS 

READ(lU)(EDGE(j) ,J=1,JB) 

IF ( ISW.NE . 1 ) WRITE( 6 ,5 ) (EDGE( J) , J=1 , JB) 

READ(IU)(FX(J),J=1, NBANDS) 

IF ( ISW.NE . 1 ) WRITE ( 6 ,6)(FX(j) ,J=1, NBANDS) 

READ(lU)(WTX(j) , J=1 , NBANDS) 

IF( ISW.NE . 1 ) WRITE ( 6,7) (WTX( j) ,  J=1 , NBANDS) 

READ(lU)(H(J) , J=1 , 64) 

I F( ISW.NE . 1 ) WRITE ( 6 , 3) 
ifCisw.eq.Dgo  TO  200 
N=NFILT/2+l 
DO  100  1=1, N, 5 
IST=I+4 

IF(IST.GT.N)IST=N 

WRITE(6,4)l,(H(ll),II=I,IST) 

100  CONTINUE 

WRITE (6, 8) 

200  ISW=1 
RETURN 

1  FORMAT (//,'  FILTER  PARAMETERS') 

2  FORMATC 5X, '  NFILT  =  ',15,'  JTYPE  =  ',15,'  NBANDS  =  ',15) 

3  FORMAT(5X, '  FILTER  COEFFICIENTS ' ) 

4  FORMAT ( 5X,I5,5E15.8) 

5  FORMAT (5X,'  EDGE  '  ,20F6.4) 

6  FORMAT( 5X, '  FX  ' ,10F6.2) 

7  FORMATC 5X, '  WTX  ’,10F7.1) 

8  FORMATC//) 

END 


1 

5 

10 

15 

20 

25 

30 

35 

40 


o  o  o  n  o 


SUBROUTINE  CONVOL(H,X,NDIM, NFILT, JTYPE) 


C 


THIS  CONVOLUTION  IS  VALID  ONLY  FOR  NFILT  =  ODD  INTEGER 
NFILT  MAX.  SET  TO  1023  -  CAN  BE  RESET  BE  REDIMENSIONING  ARRAYS 
IF  JTYPE  =  1,  EVEN  SYMMETRY  ASSUMED  FOR  REFLECTION  AT  ENDPOINTS 
IF  JTYPE  =  2,  ODD  SYMMETRY  ASSUMED  FOR  REFLECTION  AT  ENDPOINTS 


DIMENSION  H( 5 1 2 ) ,X(NDIM) ,S(1030) ,T(1030) 

IF ( JTYPE. EQ. 1 )SGN=1 . 

IF( JTYPE. EQ.2)SGN=-l • 

IF ( JTYPE. EQ.1)A=0. 

IF(JTYPE.EQ.2)A=2. 

IF( JTYPE.LT. 1. OR. JTYPE.GT.2)WRITE( 6, 1) 

IF ( JTYPE . LT . 1 . OR . JTYPE . GT . 2 ) STOP 

L=NFILT-1 

NCOEF=NFILT/2+l 

J=NCOEF 

K=  J-l 

DO  10  1=1 , NCOEF 
S(l)=SGN*X(NCOEF-I+l)+A*X(l) 

S( NFILT- 1+1 )=X( J) 

T(NCOEF+I-l)=SGN*X(NDIM-I+l)+A*X(NDIM) 

J=J-1 

10  CONTINUE 

DO  40  1=1 , NDIM 

X(I)=.0 

DO  20  J=1,K 

X( l)=X(l)+H( J)*( SGN*S( J)+S(NFILT-J+1 ) ) 

20  CONTINUE 

X( l)=X( l)+H(NCOEF )*S( NCOEF ) 

IF(l.EQ.NDIM)GO  TO  50 
DO  30  J=1 , L 
S(J)=S(J+1) 

30  CONTINUE 

if( i.le .ndim-ncoef) s(nfilt)=x( i+ncoef) 

IF ( I . GT . NDIM-NCOEF) S ( NF  ELT)=T ( I+NFILT-NDIM+ 1 ) 

40  CONTINUE 
50  RETURN 

1  FORMATC  ERROR  IN  SUB  CONVOL  -  JTYPE  NOT  EQUAL  TO  1  OR  2') 
END 
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SUBROUTINE  CONVT(M, N, ISW) 

CONVERSION  OF  DATA  (IN  COUNTS)  TO  DISPLACEMENTS 

COMMON/ I / MIN ( 6 ) , NIN(6 ) ,NTS( 24) ,NCH , IDERIV , IFILT , NPT , LDEL 

*  , NLSQ , LABEL( 10) , LABL( 3,2,24), LCNT( 2,7) 

REAL  LDEL 

COMMON/LEVl/DTIME,CT(3,2) , CT1 ( 3 , 2) ,CT2( 3 , 2) , CT3( 3 , 2 ) , CVT( 3 , 6) 

*  , TSTART , ZZ( 3 , 2) 

COMMON/LEV2C/DAT( 3 , 2 , 5000 ) , DATD ( 3 , 2 , 5000) , DATDD ( 3 , 2 , 5000 ) 

LEVEL  2 , DAT , DATD , DATDD 

ISW=ISW+1 

J=N 

IF(M.EQ. 1 .AND.N.EQ.2) J=4 
WRITE(6,1)M,N,CT(M,N) ,CVT(M,J) 

DO  100  L=1 ,NPT 

DAT(M,N,L)=DAT(M,N,L)*CT(M,N) 

IF(M.E0.1)G0  TO  100 
DAT(M,N,L)=DAT(M,N,L)*CVT(M,N) 

100  CONTINUE 

CHANGE  MICRONS  TO  CENTIMETERS 

DO  150  L=1 ,NPT 
DAT(M,N,L)=DAT(M,N,L)*l.E-4 
150  CONTINUE 

IF(M.GT.l)GO  TO  300 

IF ( MOD ( ISW, 2) .NE.O)GO  TO  300 

DO  200  L=1 , NPT 

SQ1=SQRT(1.+CVT(1,2)*DAT(M,1,L)**2+CVT(1,3)*DAT(M,2,L)**2) 
SQ2=SQRT( 1 .+CVT( 1 ,5)*DAT(M, 1 ,L)**2+CVT( 1 ,6)*DAT(M,2,L)**2) 

DAT ( M , 1 ,L)=CVT( 1 , 1 )*DAT(M, 1 ,L)/SQl 
DAT(M,2,L)=CVT(1,4)*DAT(M,2,L)/SQ2 
200  CONTINUE 
300  RETURN 

1  FORMAT(2I3 , 2F10.5) 

END 
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SUBROUTINE  DATAIN(M,N,Z, IU) 

INPUT  EXPERIMENTAL  DATA 

COMMON/ l/MIN( 6) , NIN(6) ,NTS(24) ,NCH, IDEPIV, IFILT ,NPT, LDEL 

*  ,NLSQ,LABEL( 10) ,LABL( 3 , 2, 24) , LCNT( 1 , 7 ) 

REAL  LDEL 

COMMON/LEV1  /DTIME  ,CT(3,2),CT1(3,2),CT2(.3,2),CT3(3,2),  CVT(  3 ,6) 

*  /TSTART.ZZO.Z) 

DIMENSION  Z( 5000) 

LEVEL  2, Z 

READ (IU) NTS 
DO  50  1=1,24 
LABL (M,N,l)=NTS(l) 

50  CONTINUE 

WR.ITE(6 , 3)M,N,  (LABL(M,N,  I)  ,1=1,24) 

READ( IU)NPT 
READ(IU)TSTART, DTIME 
WRITE (6 , 2 )NPT , TSTART , DTIME 
READ ( IU ) ( Z ( I ) , 1= 1 , NPT ) 

RETURN 

2  FORMAT(llO , 2F10.5) 

3  FORMAT( IX, 212 , IX, 24A2) 

END 


1 

5 

10 

15 

20 

25 

30 

35 


SUBROUTINE  DIFFER 


C 

C  CALC  FIRST  AND  SECOND  DERIVATIVES 

C 

COMMON/ l/MIN( 6) ,NIN(6) ,NTS(24) ,NCH, IDERIV , IFILT , NPT , LDEL 

*  , NLSQ , LABEL( 10  )  , LABL( 3,2,24), LCNT( 2,7) 

REAL  LDEL 

COMMON/LEVl /DTIME ,CT(3,2),CT1(3,2),CT2(3,2),CT3(3,2), CVT( 3 ,6) 

*  , TSTART , ZZ( 3 , 2) 

COMMON/LEV2C/DAT( 3 , 2 , 5000) ,DATD( 3 , 2 , 5000) , DATDD( 3 , 2 , 5000) 
LEVEL  2 , DAT , DATD , DATDD 
DO  500  1=1, NCH 
M=MIN(I) 

N=NIN(I) 

DO  100  J= 1 , NPT 

DATD ( M , N , J ) =DAT ( M , N , J ) 

100  CONTINUE 
ITYPE=2 

CALL  FILT ( DATD, M,N, NPT, ITYPF) 

DO  150  J=  1 ,  NPT 

DATD(M,N,J)=DATD(M,N,J) /DTIME 
150  CONTINUE 

WRITE (6, 1)(LABL(M,N, J) , J=1 ,24) 

ITYPE=1 

CALL  FILT ( DATD, M,N, NPT, ITYPE) 

WRITE(6, 2) (LABL(M,N, J) , J=1 , 24) 

DO  200  J=  1 ,  NPT 
DATDD(M,N,J)=DATD(M,N,J) 

200  CONTINUE 
ITYPE=2 

CALL  FILT ( DATDD, M,N, NPT, ITYPE) 

DO  250  J=1 , NPT 

DATDD (M , N , J ) =DATDD ( M , N , J ) / DTIME 
250  CONTINUE 

WRITE (6 , 1 ) (LABL(M,N, J) , J=1 , 24) 

ITYPE=1 

CALL  FILT ( DATDD, M,N, NPT, ITYPE) 

WRITE (6 , 2) (LABL(M,N, J) , J=1 , 24) 

500  CONTINUE 


RETURN 

1  FORMATC  CHANNEL  1 , 24A2 , '  HAS  BEEN  DIFFERENTIATED') 

2  FORMATC  CHANNEL  ' , 24A2 , '  HAS  BEEN  FILTERED') 

END 
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SUBROUTINE  DI S  S4 ( X , Y1 , Y2 , Y3 , NPT , LABEL ,  X4  ,  Y4  ,  N4  ) 


C 


GENERATES  PLOT  USING 


DISSPLA  (VERSION  8.2) 


DIMENSION  LABEL ( 10) ,X(NPT) ,Yl(NPT) ,Y2(NPT) , Y3(NPT) ,X4(N4) , Y4(N4) 
XORIG=X( 1) 

XSTP=' SCALE' 

XMAX=X(NPT) 

YORIG=1.E100 
YSTP=' SCALE' 

YMAX=-1 .E100 
DO  200  1=1, NPT 
IF(Y1(I) .LT.YMAX)GO  TO  100 
YMAX=Y1 ( I) 

100  IF(Y1(I).GT.Y0RIG)G0  TO  200 
YORIG=Yl (i) 

200  CONTINUE 

IF ( YORIG . EQ . YMAX ) YMAX=YMAX+ . 0 1 
WRITE( 6,1) NPT , XORIG , XMAX , YORIG , YMAX 
1  FORMAT (I5,4E15.5) 

CALL  HEIGHT (.28) 

J=J+1 

CALL  BGNPL(J) 

CALL  PHYSOR( 1 . , 1 . ) 

CALL  XINTAX 

CALL  SETCLR( 'BLACK' ) 

CALL  TITLE(LABEL(7) ,-11,'L' , 1 , '  ',5, 9. ,6. 5) 

CALL  GRAF ( XOR IG , XSTP , XMAX , YOP IG , YSTP , YMAX ) 

IMARK=0 

CALL  SETCLR( ' BLUE ' ) 

CALL  CURVE(X, Yl ,NPT , IMAFK) 

CALL  SETCLR('RED') 

CALL  DASH 

IF ( Y2 ( 1 ) . NE . 1 . E 100 ) CALL  CURVE ( X , Y2 , NPT , IMARK ) 

CALL  RE SET ( 'DASH' ) 

CALL  SETCLR( 'GREEN' ) 

CALL  DOT 

IF(Y3(l).NE.l .E 100) CALL  CURVE (X,Y3, NPT, IMARK) 

CALL  RESET('DOT') 


IMARK=-1  *•  J 

CALL  SETCLR(’RED’) 

IF(N4.GT.O) CALL  CURVE (X4,Y4,N4,IMARK) 
CALL  ENDPL(j) 

RETURN 

END 
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SUBROUTINE  FILT(Z ,M, N, NPT , ITYPE) 


C 

C 

C 

C 

C 


100 


120 


150 


200 


300 


INPUTS  FIR  COEFFICIENTS  AND  CONVOLVES  THE  DATA 

ITYPE  =  1  -  COEFFICIENTS  ARE  FOR  A  LOWPASS  OR  A  HIGHPASS  FILTER 
2  -  COEFFICIENTS  ARE  FOR  A  DIFFERENTIATOR 

DIMENSION  Z( 3, 2, 5000) ,X(5000) ,H(512) 

LEVEL  2,Z 

PI2=6. 283185307179865 
DO  100  1=1, NPT 
X(l)=Z(M,N,l) 

CONTINUE 

IF( ITYPE. E0.1)G0  TO  150 

JTYP=2 

IU=4 

CALL  COEF2(NFILT, H, IU, JTYP) 

DO  120  1=1 , NFILT 

H(I)=H(I)*PI2 

CONTINUE 

REWIND  IU 

GO  TO  200 

JTYP= 1 

IU-17 

CALL  COEF1 (NFILT, H, IU) 

REWIND  IU 

CALL  CONVOL(H,X, NPT, NFILT, JTYP) 

DO  300  1=1, NPT 
Z(M,N,I)=X(I) 

CONTINUE 

RETURN 

END 
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SUBROUTINE  ORIENT 


C 

C 


C 

C 


READ  CONTROL  PARAMETERS 
READ  DATA  INTO  APPROPRIATE  ARRAYS 
CONVERT  DATA,  IF  DESIRED 
FILTER  DATA,  IF  DESIRED 

COMMON/ l/MIN( 6) ,NIN(6) , NTS (24) ,NCH, IDE R IV,  IFILT,NPT,LDEL 

*  , NLSQ , LABEL( 10 ) , LABL( 3,2,24), LCNT( 2,7) 

REAL  LDEL 

COMMON/LEVl /DTIME ,CT(3,2) , CT1 ( 3 , 2 ) ,CT2(3,2) ,CT3(3,2) , CVT( 3 ,6) 

*  , TSTART , ZZ( 3 , 2) 

C0MM0N/LEV2C/DAT( 3,2, 5000) ,DATD(3 , 2 , 5000) ,DATDD( 3 , 2 , 5000) 
LEVEL  2 , DAT , DATD , DATDD 
COMMON/LEV2D/Z( 5000) 

LEVEL  2,Z 

READ( 5 , 1 )nch, icvt , ifilt,  ideriv 

WRITE ( 6 , 6 ) NCH , ICVT , IFILT , IDERI V 
IF( IDERIV. EQ.0)JL=3 
IF( IDERIV. EQ.1)JL=1 
WRITE (6, 5) 

DO  150  J=  1 ,  JL 

DO  100  1=1, NCH 

READ(5 , 1)MIN( I) ,NIN(l) , IU 

WRITE ( 6 , 1 ) MIN ( I ) , NIN( I ) , IU 

M=MIN( I) 

N=NIN(I) 

CALL  DATAIN ( M , N , Z , IU ) 

IF(J.EQ.2)G0  TO  70 
IF(J.EQ.3)GO  TO  90 
DO  60  K=  1 ,  NPT 
DAT(M,N,K)=Z(K) 

60  CONTINUE 
GO  TO  100 

70  DO  80  K=1 , NPT 

DATD (M,N,K)=Z(K) 

80  CONTINUE 
GO  TO  100 

90  DO  R5  K= 1 , NPT 
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DATDD(M,N,K)=Z(K) 

95  CONTINUE 
100  CONTINUE 
150  CONTINUE 

IF(ICVT.EO.O)GO  TO  325 
WRITE (6, 7) 

READ (5 , 3) (CVT( 1 , J) , J=1 ,6) 

WRITE (6, 3) (CVT( 1 , J) , J=1 ,6) 

READ(5,3)(CVT(2,J) , J=1 ,2) 

WRITE (6 , 3) (CVT( 2 , J) , J=1 , 2) 

RE AD (5 ,  3)  ! CVT( 3 , J) , J=1 , 2) 

WRITE(6,3)(CVT(3 , J) , J=1 ,2) 

DO  300  1-1,3 
DO  300  J=1 , 2 

READ (5 ,2)M,N,CTl!M,N) ,CT2(M,N) ,CT3(M,N) 
CT(M,N)=CT1(M,N)/CT2(M,N)/CT3(M,N) 

WRITE (6 , 2)M,N,CTl!M,N) ,CT2(M,N) ,CT3(M,N) ,CT(M,N) 

300  CONTINUE 
325  ISW=0 

DO  400  1=1, NCH 
M=MIN( I ) 

N=NIN(I) 

IF ( ICVT . EQ . 0 ) GO  TO  350 
CALL  CONVT (M,N,ISW) 

350  IF( IFILT .EQ.O)GO  TO  400 
'  ITYPE=1 

CALL  FILT(DAT,M,N,NPT, ITYPE) 

WRITE (6, 4) (LABL(M,N, J) , J=1 ,3) 

400  CONTINUE 
RETURN 

1  F0RMAT(4I3) 

2  FORMAT ( 213 ,4F10. 5) 

3  FORMAT! 8F10. 5) 

4  FORMAT! 1  CHANNEL  ’,3A2,’  HAS  BEEN  FILTERED') 

5  FORMAT!//,'  INPUT  DATA  CHANNELS') 

6  FORMAT!'  NCH  =  ’,13,'  ICVT  =  ',13,'  IFILT  =  ',13,'  IDERIV  =  ',13) 

7  FORMAT!//,'  CALIBRATION  CONSTANTS') 

END 
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SUBROUTINE  PLOT4(DAT) 


PUTS  3-DIMENSIONAL  ARRAY  INTO  1-DIMENSIONAL  ARRAYS  FOR  PLOTTING 

DIMENSION  DAT (3 , 2 , 5000) 

LEVEL  2, DAT 

COMMON/I/MINC6) ,NIN(6) ,NTS ( 24) , NCH , IDERIV , IFILT , NPT , LDEL 

*  ,NLSQ,LABEL(10) ,LABL( 3 , 2 , 24) , LCNT(2,7) 

REAL  LDEL 

COMMON  DUM(l) , 

*  X(5000) ,Y1 ( 5000) ,Y2(5000) ,Y3(5000) ,Y4(5000) ,Y5(5000) ,Y6(5000) 
DUM(1)=1.E100 

NPL1’=2500 
.  N4=0 

ENCODE ( 11,1, LABEL ( 7 ) ) (LABL( 1 , l , J) , J=1 , 3) 

DO  100  J=1 , NPLT 
X(J)=J 

Y1(J)=DAT(MIN(1),NIN(1),J) 

IF(MIN( 3) .NE.0)Y2(J)=DAT(MIN(3) ,NIN(3) , J) 

IF(MIN(5) .NE.0)Y3(J)=DAT(MIN(5) ,NIN(5) ,J) 

Y4(J)=DAT(MIN(2) ,NIN(2) , J) 

IF(MIN(4) ,NE.0)Y5(J)=DAT(MIN(4) ,NIN(4) , J) 
IF(MIN(6).NE.0)Y6(J)=DAT(MIN(6),NIN(6),J) 

100  CONTINUE 

IF(MIN(3) .E0.0)Y2(l)=DUM(l) 

IF(MIN(4) ,E0.0)Y5(1)=DUM(1) 

IF(MIN(5).EO.O)Y3(l)=DUM(l) 

IF(MIN(6).EQ.0)Y6(1)=DUM(1) 

ENCODE( 11 , 1 , LABEL (7 ) ) (LABL( 1 , 1 , J) , J=1 ,3) 

CALL  DISS4(X,Y1,Y2,Y3, NPLT, LABEL, X4,Y4,N4) 

ENCODE ( 11,1 , LABEL ( 7 ) ) (LABL( 1 , 2 , J) , J=1 , 3) 

CALL  DISS4(X,Y4,Y5,Y6, NPLT , LABEL , X4 , Y4 , N4 ) 

RETURN 

1  FORMAT( 3A2 , 5X) 

END 
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SUBROUTINE  WRDAT 

WRITES  ALL  PROCESSED  DATA  ONTO  OUTPUT  FILE  TO  BE  SAVED  FOR  LATER  USE 

COMMON/LEV2C/DAT( 3 , 2 , 5000) ,DATD( 3,2, 5000) ,DATDD( 3 , 2 , 5000) 

LEVEL  2 , DAT , DATD , DATDD 

COMMON /l/MIN(6), NIN ( 6 ) ,NTS(24) ,NCH, IDERIV, IFILT.NPT.LDEL 

*  , NLSQ,LABEL( 10) ,LABL( 3 , 2 , 24) ,LCNT(2,7) 

REAL  LDEL 

COMMON/LEVI /DTIME ,CT( 3 , 2) ,CT1(3,2) ,CT2(3,2) ,CT3(3,2) ,CVT(3,6) 

*  , TSTART , ZZ ( 3 , 2 ) 

IF(IFILT.EQ.0)WRITE(6,2) 

IF( IFILT.NE.O) WRITE( 6 , 5) 

DO  200  1=1, NCH 

DO  100  J=1 , 24 

NTS( J)=LABL(MIN( I) ,NIN( I) , J) 

100  CONTINUE 

WRITE( 18)NTS 
WRITE( 6,1) NTS 
WRITE( 18 ) NPT 
WRITE(18)TSTART, DTIME 

WRITE ( 18) ( DAT(MIN( I) , NIN( I) , J) , J=1 ,NPT) 

200  CONTINUE 

IF( IDERIV. EQ.0)GO  TO  700 
WRITE (6, 3) 

DO  400  1=1, NCH 
DO  300  J=1 , 24 

NTS(J)=LABL(MIN(I),NIN(I),J) 

300  CONTINUE 

WRITE( 18)NTS 

WRITE (6,1) NTS 

WRITE( 18)NPT 

WRITE ( 1 8 ) TSTART , DTIME 

WRITE ( 18) (DATD(MIN( I) ,NIN( I) , J) , J=1 ,NPT) 

400  CONTINUE 

WRITE(6 ,4) 
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DO  600  1=1, NCH 
DO  500  J=1 , 24 

NTS(J)=LABL(MIN(I) ,NIN(I) ,J) 

500  CONTINUE 

WRITE( 18)NTS 

WRITE (6,1) NTS 

WRITE(18)NPT 

WRITE ( 1 8 ) TSTART , DTIME 

WRITE( 18)(DATDD(MIN(I) ,NIN(I) , J) , J=1 ,NPT) 
600  CONTINUE 
700  RETURN 

1  FORMAT(5X, '  WRITING  ON  FILE  -  ' ,24A2) 

2  FORMAT(//,'  INPUT  DATA  -  UNFILTERED’) 

3  FORMATC  FIRST  DERIVATIVE  -  FILTERED') 

4  FORMATC  SECOND  DERIVATIVE  -  FILTERED') 

5  FORMAT(//,'  INPUT  DATA  -  FILTERED') 

END 
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APPENDIX  C 


C.1.3.  Card  Image  Formats  That  Are  Required  For  The  Input  Data  Section  Of  MFA  File  MUZMO40/UN-BOOTS 

TABLE  C.l.  INPUT  DATA  FORMATS  FOR  TAPE  UNIT  5 


C71 


Card 

Condition* 

Column 

Format 

Variable 

Calling  Subroutine 

Description 

i 

1-3 

13 

IPLTl 

MAIN 

If  IPLT1=1,  displacements  are  plotted 

4-6 

13 

IPLT2 

If  IPLT2=1,  velocities  are  plotted 

7-9 

13 

IPLT3 

If  IPLT3=1,  accelerations  are  plotted 

10-12 

13 

IPRT 

If  IPRT=1 ,  displacements,  velocities 
and  accelerations  are  written  on  a  file 

on  tape  unit  18. 

2 

1-3 

13 

NCH 

ORIENT 

Number  of  data  sets  to  be  inputted 

4-6 

13 

ICVT 

If  ICVT  ^0,  data  conversion  is  done 

7-9 

13 

IFILT 

If  IFILT  ^0,  displacements  are  filtered 

10-12 

13 

IDERIV 

If  IDERIV  ^0,  differentiation  is  done 

3 

1-3 

13 

MIN(I) 

ORIENT 

First  index  for  placement  of  data  into 
array  DAT 

4-6 

13 

MIN(I) 

Second  index  for  placement  of  data  into 
array  DAT 

7-9 

13 

IU 

Tape  unit  number  on  which  data  file 
resides 

4 

Repeat  card 

3  until  NCH 

cards  are 

inputted . 

5.  Repeat  sets  of  cards  3  and  4  for  JL  sets:  JL-3  if  IFILT  -  0  or  JL  1  if  IFILT  1* 


*Omit  the  card  if  condition  is  not  met. 


TABLE  C.l.  INPUT  DATA  FORMATS  FOR  TAPE  UNIT  5  (continued) 


cn 

GO 


Card 

Condition 

Column 

Format 

6 

ICVT  1 

0 

1-10 

F10.5 

11-10 

F10.5 

CVT( 1,2) 

21-20 

F10.5 

CVT( 1,3) 

31-40 

F10.5 

CVT( 1,4) 

41-30 

F10.3 

CVT( 1,5) 

51-60 

F10.5 

CVT( 1,6) 

7 

ICVT  f 

0 

1-10 

F10.5 

11-20 

F10.5 

CVT(2 ,2) 

8 

ICVT  i 

0 

1-10 

F10.5 

11-20 

F10.5 

CVT( 3,2) 

9** 

ICVT  1 

0 

1-3 

13 

Variable  Calling  Subroutine 

CVT( 1,1)  ORIENT 


CVT( 2,1)  ORIENT 

CVT( 3,1)  ORIENT 

M  ORIENT 


Description 

Coefficients  for  conversion 
of  data  of  sensor  1 
[conversion  is  nonlinear, 
as  in  Eq .  (3.2)] 

Coefficients  for  conversion 
of  data  of  sensor  2 
[conversion  is  linear] 

Coefficients  for  conversion 
of  data  of  sensor  3 
[conversion  is  linear] 

First  index  of  arrays  CT1,CT2,CT3 
Second  index  of  arrays  CT1,CT2,CT3 

Transducer  calibration  factors 


4-6 

7-16 

17-26 

27-36 

10  ICVT  f  0 


13  N 

F10.5  CT1 (M, N) 

F10.5  CT2(M,N) 

F10.5  CT3(M,N) 

Repeat  card  9  until 


6  cards  are  inputted 


**Since  these  data  sets  had  not  been  previously  converted  into  voltages  so  this  input  card  was  incorporated 
into  the  calculations  such  that  CT(M,N)=CT1(M,N)/CT2(M,N)/CT3(M>N) .  The  data  are  then  multiplied  by 
the  appropriate  CT  before  the  other  conversion  is  done* 


C.1.4.  Output 


The  normal  printed  output  of  this  program  is  mainly  informative  and  diagnostic.  The  output  consists  of 

1.  Input  variables; 

2.  Identification  of  input  data  channels; 

3.  Calibration  constants; 

4.  Lowpass  filter  coefficients  (when  used  first  time); 

5.  Differentiator  coefficients  (when  used  first  time);  and 

6.  Path  each  channel  follows  as  it  is  differentiated,  filtered,  and/or  written  on  an  output  file. 

The  plotted  output  consists  of  sets  of  two  plots:  one  for  horizontal  data  and  one  for  vertical.  The  number 
of  curves  on  each  plot  is  determined  by  the  array  MIN.  The  number  of  sets  plotted  is  determined  by  the 
variables  IPLT1,  IPLT2 ,  and  IPLT3 • 


C.2.  Program  MUZPRED 

C.2.1.  Listing  Of  Job  Control  Language  Of  MFA  File  MOZPRED/DH=BOOTS 

BOOTS, T200. 

USER , BOOTS , XXXXXXX . 

CHARGE, PDXXX, PD. 

ATTACH, DISSPL9/UN=DISSPLA. 

ATTACH, IMSL/UN=LIBPARY. 

LIBRARY, DISSPL9,IMSL. 

FTN,L=0. 

GET, A=DIFFATR. 

FTN , I=A,L=0 . 

ATTACH , TAPE 1 =DVA05 . 

GET,TAPE2=DIFF1. 

GET , TAPE3=DISPL40 . 

ATTACH , TAPE4=DVA0  5U . 

GET, TAPE 7=FECL40. 

GET , TAPE8=LPF67 . 

GET , TAPE 1 0=UVECT40 . 

LGO. 

REPLACE, TAPE9=MUZPT05. 

REPLACE, META=DISS05R. 

EXIT. 

REPLACE ,META=DISSERR . 


C.2.2.  Listing  Of  FORTRAN  Program 
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5 


C 

C 

C 

C 

C 
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PROGRAM  DISS9 ( INPUT , OUTPUT , TAPE5= INPUT , TAPF6=OUTPUT , TAPE 1 , TAPE2 , 

*  TAPE 3 , TAPE4 , TAPE 7 , TAPE 8 , TAPE 9 , TAPE  1 0 , META , TAPE 1 3 ) 

ANALYSIS  OF  MUZZLE  MOTION  DATA  -  PROGRAM  2 

PLOTTING  IN  COLOR  USING  DISSPLA  9.0  ON  A  TEKTRONIX  4107  TERMINAL 

COMMON/ LEVI /DAT (6 ,2,2500) , LABL( 3,2,24) ,MIN( 9) ,NIN( 9 ) ,TSTART,DTIME , 

*  NPT, NPLT,LKEY 
DIMENSION  LABEL ( 10) 

LABEL( 1 )=10H  BOOTS 
LABF,L(2)  =  10HBLDG.390 
LABEL( 3)=10HX  6121 
LABEL(4)=10HDISSPLA9.0 
CALL  COMPRS 

CALL  SETDEV( 13 , 13) 

WRITE ( 6,1) 

RE AD ( 5 , 2 ) NCH , NPLOTP , NPLOTC , NPLOTU , IPR I NT 
IF(EOF(5).NE..O)GO  TO  200 

WRITE ( 6 , 7 ) NCH , NPLOTP , NPLOTC , NPLOTU , IPR I NT 

DO  110  1=1, NCH 

READ( 5 , 2)MIN( I) ,NIN( I) , IU 

WRITE ( 6 , 6)MIN( I) , NIN( I) , IU 

IF( IU.EQ.0) IU=1 

CALL  DATAIN(MIN(I),NIN(I),IU) 

110  CONTINUE 

READ(5,3)LKEY 
WR ITE (6,4) LKEY 
CALL  PREDCT( NPLOTP) 

CALL  CALCZ ( NPLOTC , NPLOTU , IPR INT ) 

200  CALL  DONE PL 
STOP 

1  FORMAT(lHl) 

2  FORMAT( 2013) 

3  FORMAT( 1615 ) 

4  FORMAT (//,'  LKEY  =  ’,15) 

6  FORMAT('  MIN,  NIN,  IU  -’,2013) 

7  FORMATC  NCH,  NPLOTP,  NPLOTC,  NPLOTU,  IPRINT  -’,2013) 

END 
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SUBROUTINE  CALCZ(NPLOTC , NPLOTU, IPRINT) 

C  ' 

C  THIS-  SUBROUTINE  SCALES  A  SET  OF  RECOIL  MOTION  DATA  FROM  A  SIMILAR 

C  FIRING  TO  THE  TIME  OF  THIS  FIRING,  CALCULATES  THE  DISPLACEMENT  IN 

C  THE  Z-DIRECTION  OF  EACH  COMPONENT  (HORIZONTAL  AND  VERTICAL)  AT  EACH 

C  DISCRETE  TIME,  DIFFERENTIATES  THIS  RESULTING  DISPLACEMENT,  AND 

C  PLOTS  THE  RESULTS  IF  NPLOTC  .NE.  0. 

C 

DIMENSION  C(10) ,MDAT(9) ,NDAT(9) ,AA(3,3) ,BB( 3 ) , COEF( 3 ) , IWK(50) , 

*  WRK(50) ,Z(3,2) ,ZZ(2500) 

COMMON/LEVI /DAT (6, 2, 2 5 00) ,LABL( 3 , 2 , 24) ,MIN(9) ,NIN(9) , TSTART,DTIME, 

*  NPT,NPLT,LKEY 
C 

LZM  =  MODEL  L( START  OF  MOTION) 

LPM  =  MODEL  L( PEAK  PRESSURE) 

LEM  =  MODEL  L(SHOT  EJECTION) 

LPD  =  DATA  L(PEAK  PRESSURE) 

LED  =  DATA  L(SHOT  EJECTION) 

LZD  IS  THEN  DETERMINED  BY  SCALING  THE  TWO  SETS  OF  PARAMETERS 
C 

READ(5, 1)LZM,LPM,LEM 
READ( 5 , 1 )LPD, LED 

RAT=FLOAT ( LEM-LZM) /FLOAT ( LEM-LPM) 

LZD=  FLO AT ( LED ) -RAT*  FLOAT ( LED-LPD ) 

RLSPAN=LED-LZD 
WRITE (6, 6) 

WRITE (6,1) LZM , LPM, LEM 
WRITE ( 6 , 1 ) LZD , LPD , LED 
READ( 5 , 1 )MPOLY 
MPP=MPOLY+l 
WRITE( 6,7) 

READ ( 5 ,2)(C(M) ,M=1 ,MPP) 

WRITE(6,2)(C(M),M=1,MPP) 

RVMAX= . 0 
DO  200  L=  1 ,  NPT 
DAT( 6 , 1 ,L)=.0 
DAT(6,2,L)=.0 
FL=FLOAT(L-LZD ) / RLSPAN 
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IF(L.LE.LZD)GO  TO  200 
DO  100  M=1,MPP 

DAT(6 , 1 ,L)=DAT(6 , 1 ,L)+C(M)*FL**(M-1 ) 

DAT( 6 , 2 ,L)=DAT( 6 j  2 ,L)+C(M)*FL**M/FLOAT(M) 

100  CONTINUE 

DAT(6,2,L)=DAT(6,2,L)*RLSPAN*DTIME 
IF(DAT( 6 , 1 ,L) . LT.RVMAX)GO  TO  260 
RVMAX=DAT( 6 , 1 , L) 

200  CONTINUE 

260  DO  270  LL=L,NPT 

DAT(6, 1 ,LL)=RVMAX 

DAT  ( 6 , 2 , LL )  =DAT  ( 6 , 2 , LL- 1 )  +SVMAX*DT  IME 
270  CONTINUE 

READ(5,2)Z(1,1),Z(1,2),Z(2,1),Z(3,1) 

Z(2,2)=Z(2,1) 

Z(3,2)=Z(3,1) 

WRITE (6 ,8) 

WRITE (6, 2) Z 
C 

C  TEMPORARY  -  THE  VARIABLE  CONV  IS  NECESSARY  TO  CHANGE  THE  RECOIL 

C  DATA  FROM  MILLIMETERS  TO  CENTIMETERS 

C 

CONV=.l 

C 

DO  450  L=1,NPT 
DO  350  K=1 , 3 
AA(X,1)=1. 

AA(K,2)=Z(K,l)+CONV*DAT(6,2,L) 

AA(K,3)=AA(K,2)**4 
bb(k)=dat(k, 1 ,l) 

350  CONTINUE 
TOL=.0 
KBASIS=0 

CALL  LLSQF( AA, 3 , 3 , 3 , BB , TOL , KBASIS , COEF , WRK , IWK , IER) 

DAT(4, 1 ,L)=C0EF( 1) 

DAT(5,l,L)=COEF(2) 

DAT( 1,1 ,L)=COEF(3) 

DO  400  K=1 ,3 
AA(K j 1 )=1 . 
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AA(K,2)=Z(K,?)+CONV*DAT(6,2,L) 

AA(E,3)=AA(K,2)**4 

BB(K)=DAT(K,2,L) 

400  CONTINUE 
T0L=.0 
KBASIS=0 

CALL  LLSOE( AA , 3 , 3 , 3 , BB , TOL , ERAS 18 , COEF , WRK , IWK , IFF) 
DAT(4,2JL)=COEF(l) 

DAT(5,2,L)=COEFC2) 

DAT(l,2,L)=COEF(3) 

450  CONTINUE 

READ(5,4)NZREF 
WRITE (6 , 10)NZREF 
DO  480  1=1 , NZREF 
READ ( 5 , 2)ZREF 
WR ITE( 6 ,  11 ) ZREF 
DO  454  L= 1 , NPT 

DAT (2 , 1 ,L)=DAT(4 , 1 ,L)+OAT(5 , 1 ,L)*ZFEF+DAT( 1 , 1 ,L)*ZFEF**4 
DAT( 2 , 2 ,L)=OAT(4 , 2 ,L)+OAT( 5 , 2 ,L)*ZREF+DAT( 1 , 2 ,L)*ZFEF**4 
DAT ( 3 , 1  ,L)=1'ATC  5 , 1  ,L)+4.*DAT(  1 , 1  ,L)*ZREF**3 
DAT( 3 , 2 ,L)=DAT( 5 , 2  jL)+4.*DAT( l , 2 ,L)*ZREF**3 

454  CONTINUE 

TKE Y=DT IME*FLOAT ( LKE Y ) 

TST=-1 .+TKEY 
TSP=  2.+TKEY 

LST=(TST-TSTART)/DTIMF,+  1  .^-lO+l . 

LSP=(TSP-TSTART)/ DTIME+1 .E-10 
IU=9 

IF( IPRINT.EO.O)GO  TO  455 
NARR=3 

READ ( 5 , 4 ) ( MDAT ( II) ,NDAT( II) , 11=] ,NARR) 

WRITE ( 6 , 9) ( MDATf II) ,NDAT( II) , 11=1 , NARR) 

CALL  WRFILE( DAT, MDAT , NDAT ,NARR , TST , LST , LSP , LABI, , DTIME , IU) 

455  DO  456  L=1,NPT 
zz(l)=dat( 1 , 1 ,l) 

456  CONTINUE 

CALL  UNITV(NPLOTU,  IPRINT,TST,LST,LSP,  III) 

DO  457  L= 1 , N^T 

dat(i,i,l)=zz(l) 
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457  CONTINUE 


C 

C  IUD  =  INPUT  TAPE  UNIT  FOR  DIFFERENTIATOR 

C  IUF  =  INPUT  TAPE  UNIT  FOR  FIR  FILTER 

C  IF  IUD  OR  IUF  =  0,  THAT  OPERATION  IS  NOT  DONE. 

C 

IUD=2 

IUF=8 

DO  460  L=1 , NPT 
ZZ(L)=DAT( 2 , 1 , L) 

460  CONTINUE 

CALL  DIFFER( ZZ , NPT , DTIME , IUD , IUF) 

DO  462  L=1 , NPT 
DAT(3,1,L)=ZZ(L) 

ZZ(L)=DAT(2,2,L) 

462  CONTINUE 

CALL  DIFFER ( ZZ , NPT , DTIME , IUD , IUF ) 

DO  464  L=1 , NPT 
DAT (3 , 2 ,L)=ZZ(L) 

464  CONTINUE 

IF ( NPLOTC . EQ . 0 ) GO  TO  480 
WRITE( 6 ,3) 

NCURV=1 
IUPLT=7 
REWIND  IUPLT 

READ( IUPLT , 5 ) TIMEL , TIMER 
READ( iuplt,4)laxis, icolor 
DO  470  11=1 .NPLOTC 
ITYPE=0 

IF( II .EQ. 3. OR. II. EO. 6. OR. II. EQ. 13) ITYPE=1 
N4=l 

IF( ITYPE.EQ. 1 )N4=2 
NCRV=NCURV*( ITYPE+1 ) 

READ ( IUPLT , 4 ) ( MDAT ( K ) , NDAT ( K ) ,K=1 ,NCRV) 
WRITE(6,9)(MDAT(K) ,NDAT(K) ,K=I,NCRV) 

ISC=0 

IF( II.EQ.4)ISC=1 
IF( ITYPE.EQ. 1 ) ISC=1 
IF(II.EQ.6)TIMEL=-1. 

IF( II .EQ.6 )TIMER=  1. 


160 


165 


170 


CALI,  DATA4 ( NCURV ,MDAT , NDAT ,  ISC , ITYPE , LAXIS , ICOLOR , N4 , TIMEL .TIMER , 
*  IUPLT) 

470  CONTINUE 
480  CONTINUE 
500  RETURN 

1  FORMAT( 1615 ) 

2  FORMAT(8E10.3) 

3  FORMAT (//,'  PLOTTING  IN  SUBROUTINE  CALCZ') 

4  FORMAT( 2013) 

5  F0RMAT(2F5.0) 

6  FORMAT (// , '  SUBROUTINE  CALCZ  -  SCALING  PARAMETERS') 

7  FORMATC  COEFFICIENTS  OF  RECOIL  MOTION  CURVE') 

8  FORMATC  Z-AFRAY') 

9  FORMATC/, '  "MDAT,  NDAT"  PAIRS  -',2013) 

10  FORMATC  NZREF  =',I3) 

11  FORMATC/'  ZREF  =  ' ,F10.5) 

END 


<T> 
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SUBROUTINE  DATAIN(M,N, IU) 


C 

C 

C 


50 


C 

C 

C 

C 


100 

1 

2 

3 


INPUT  EXPERIMENTAL  DATA 

COMMON/LEVI /DAT ( 6 ,2,2500), LABL( 3 ,2,24) ,MIN(9) ,NIN(9) ,TSTART, DTIME, 

NPT ,  NPLT ,  I, KEY 

DIMENSION  NTS( 24) , Z(4100) 

READ(IU)NTS 
DO  50  1=1,24 
LABL(M,N,I)=NTS(I) 

CONTINUE 

WRITE(6,3)M,N,(LABL(M,N,I) ,1=1,24) 

READ ( IU  ) NDATA 
READ( IU)TSTART, DTIME 
WRITE(6, 1)NDATA,TSTART,DTIME 
READ( IU)(Z(l),I=l, NDATA) 

NPT= NDATA 

NECESSARY  BECAUSE  OF  CENTRAL  MEMORY  SIZE  RESTRICTIONS 

IF(NPT.GT.2500)NPT=2500 

NPLT=NPT 

DO  100  1=1, NPT 

DAT(M,N,I)=Z(I) 

CONTINUE 

RETURN 

FORMAT ( '  NDATA,  TSTART,  DTIME  ,I10,2F10.5) 

FORMAT (I10,2F10.5) 

FORMAT ( IX, 212, IX, 24 A2) 

END 
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SUBROUTINE  DATA4 C NCURV , MDAT , NDAT , ISC , ITYPE , LAXIS , ICOLOR , N4 , TIMEL , 

*  TIMER, IUPLT) 

C 

THIS  SUBROUTINE  PREPARES  THE  DATA  ARRAYS  IN  PROPER  FASHION  SO 
THAT  A  WIDE  VARIETY  OF  PLOTS  CAN  BE  GENERATED  USING  THE  SAME 
PLOTTING  SUBROUTINE. 

DIMENSION  MDAT ( 9 ) ,NDAT(9) ,X( 2500) ,Y1 (2500) ,Y2(2500) ,Y3(2500) , 

*  DUM( 1),X4(8) ,Y4(8) ,LXNAME(2) , LYNAME(2) 

COMMON/LEVl/DAT(6 , 2,2500),LABL(3,2,24) ,MIN( 9) , NIN(9) ,TSTART, DTIME , 

*  NPT,NPLT,LKEY 
X4  ( 1 )  = . 0 

IF(N4.EQ.O)GO  TO  40 
RE AD ( IUPLT, 3) ( (X4( I) , Y4( I) ) , 1=1 ,N4) 

DO  30  1=1, N4 
X4( l)=X4( l)*DTIME 
30  CONTINUE 

WRITE (6 , 6) ( (X4(l) , Y4( I) ) , 1=1 ,N4) 

40  DUM( 1 )=1 .E100 
TST=X4( 1 )+TIMEL 
TSP=X4 ( 1 ) +TIMER 

LST=(TST-TSTART)/DTIME+1 . E-10+1 . 

LSP=(TSP-TSTART)/DTIME+1.E-10 
NPLOT=LSP-LST+ 1 
TSTRT=-X4(1) 

X4  ( 1 )  = .  0 


SPECIAL  CONDITION  -  IF 
C  COORDINATES 

C 

IF(N4.EQ.2)N4=0 


C 


J=0 


N4 


2, 


DO  100  L=LST,LSP 
J=J+1 

IF( ITYPE. EQ.l)GO  TO  50 
X(J)=TSTRT+DTIME*FL0AT(L-1) 
Y1 ( J)=DAT(MDAT( 1 ) ,NDAT( 1 ) ,L) 


THE  PLOT  ORIGIN  IS  MADE  AT  THESE 
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IF(NCURV.GE.2)Y2(J)=DAT(MDAT(2),NDAT(2) ,L) 

IF  C  NCUR V . GE . 3 ) Y3  C  J ) =DAT ( MDAT ( 3 ) , NDAT (  3  )  ,  L  ) 

GO  TO  100 

50  X ( J ) =DAT ( MDAT ( 1 )  , NDAT ( 1 )  ,  L ) 

Yl ( J ) =DAT ( MDAT ( 2 ) , NDAT C  2  )  ,  L  ) 

IF(NCURV.GE.2)Y2(j)=DAT(MDAT(3) ,NDAT(3) ,L) 

IFCNCURV.GE . 2) Y3 ( J)=DAT(MDAT(4) ,NDAT(4) ,L) 

IF(N4.EQ.0)X(j)=X(j) -DAT ( MDAT ( 1 )  , NDAT ( 1 ) , LKEY ) 

IF(N4.EQ.0)Y1(J)=Y1 ( J)-DAT(MDAT( 2) , NDAT (2 ) ,LKEY) 

100  CONTINUE 

IF(N4.EQ.0)WRITE(6,5)DAT(MDAT(1) ,NDAT(1) ,LKEY) ,DAT(MDAT(2) ,NDAT(2) 

*  , LKEY ) 

IF(NCURV.LT. 2)Y2(1)=DUM(1) 

IF(NCURV.LT. 3.AND.ITYPE.EQ.O)Y3(l)=DUM(l) 

ENCODE  C 10 , 1 , TITLE) C LABL C  MDAT  Cl), NDAT  C 1 ) , J ) , J= 1 , 3 ) 

WRITE  C  6 , 2 ) LST , LSP , X ( 1 ) , X  C  NPLOT) 

LXNAMECl)='  ' 

LYNAMECl)='  ' 

LXNAMEC2)='  ' 

LYNAMEC2)='  ' 

IF  C  LAXI S . EQ . 1 ) READ  C IUPLT , 4 ) LXNAME , LYNAME 
IF  C  LAXIS . EQ . 1 ) WRITE  C  6 , 4 ) LXNAME , LYNAME 

CALL  DISS4CX,Y1,Y2,Y3, NPLOT , TITLE, X4,Y4,N4, ISC, ITYPE, LXNAME , LYNAME 

*  , ICOLOR, IUPLT) 

RETURN 

1  FORMAT C 3 A2,4X) 

2  FORMATC  PLOT  LIMITS  -  '  ,2I5,2F10.4) 

3  FORMATC 8F 10.0) 

4  FORMATC 8A10) 

5  FORMATC  COORDINATES  AT  SHOT  EJECTION  -  ' , 2 CE12 . 5 , 3X) ) 

6  FORMATC  SPECIAL  COORDINATES' ,8F10. 3) 

END 
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SUBROUTINE  DISS4(X, Y1 ,Y2 ,Y3 , NPT , TITLE , X4 ,Y4 ,N4 , ISC , ITYPE , LXNAME 
*  LYNAME , ICOLOR, IUPLT) 


C 

C 

C 

C 

C 

C 

C 

C 


ISC=0  SELF-SCALE 

=1  READ  IN  XORIG,  XMAX,  YORIG,  YMAX 


ITYPE=0  NORMAL  PLOT 
=1  POLAR  PLOT 


ICOLOR=0  BLACK/WHITE 
=1  COLOR 


DIMENSION  X(NPT) ,Y1 (NPT) , Y2(NPT) ,Y3(NPT) ,X4(N4) ,Y4(N4) 
*  ,RAT(  10) ,LXNAME(2) ,LYNAME(2) 

IF(ISC.NE.O)GO  TO  300 
XORIG=X ( 1 ) 

XSTP=' SCALE' 

XMAX=X ( NPT ) 

YORIG=1.E100 
YSTP= ' SCALE ' 

YMAX=-1.E100 
DO  200  1=1, NPT 
IF(Yl(l).LT. YMAX ) GO  TO  100 
YMAX=Yl ( I) 

100  IF(Y1 ( I) .GT. YORIG) GO  TO  200 
Y0RIG=Y1 ( I) 

200  CONTINUE 


C 

C 

C 


SPECIAL  FOR  THIS  SET  OF  DATA 


IF(YORIG.LT. .0)YORIG= YORIG* 1.1 


C 


GO  TO  400 

300  READ (IUPLT, l)XORIG, XMAX, XSTP, YORIG, YMAX, YSTP 
400  WRITE (6, 2) XORIG, XSTP, XMAX, YORIG, YSTP, YMAX 
JJ=JJ+1 

CALL  BASALFC ' STANDARD' ) 

CALL  MIXALF( ' SPECIAL' ) 

CALL  PHYSOR( 1.5,1.) 
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CALL  PAGE ( 1 1 . , 8 . 5) 

CALL  SETCLR ( 'BLACK' ) 

TF( ITYPE . E0 . 0 ) CALL  AREA2D( 8 .5,6.0) 

IF( ITYPE .EQ. 1 ) CALL  AREA2D(6.0,6.0) 

CALL  HEIGHT( . 25 ) 

CALL  HEADIN(TITLE, 10, 1.5,1) 

CALL  XNAME ( LXNAME ,13) 

CALL  YNAME ( LYNAME ,13) 

CALL  GR AF ( XOR IG , XSTP , XMAX , YOR IG , YST P , YMAX) 

IF ( ITYPE . EQ • 0 ) IMARK=0 
IF( ITYPE. E0.1)IMARK=100 
IF(ITYPE.EQ. 1 .AND.N4.EO.O)IMARK=25 
IF( ICOLOR . EO . 1 ) CALL  SETCLR ( ' BLUE ' ) 

CALL  CURVE(X,  Yl ,  NPT,  IMARK) 

IF( ITYPE. EO.O) GO  TO  600 
DO  500  11=1 ,NPT, 50 

CALL  RLVEC(X(ll),Yl(ll),X(lI+5),Yl(lT+5) ,2331) 

500  CONTINUE 

600  IF(  ICOLOR. EO.UCALL  SETCLR(  '  RED '  ) 

TL=.5 
NMRK=2 
RAT( 1 )=6 . 

RAT(2)=4. 

IF ( ICOLOR .EQ.O) CALL  MRSCOD ( TL , NMRK , RAT ) 

IF(Y2(1) .NE. 1.E100. AND. ITYPE. EO.O)CALL  CURVE (X,Y2,NPT, IMARK) 
IF ( Y2  C 1 ) . NE . 1 . E 1 00 . AND . ITYPE . EO • 1 ) CALL  CURVE ( Y2 , Y3 , NPT , IMARK ) 
CALL  RE SET ('DOT') 

IF ( ICOLOR . EO . 1 ) CALL  SETCLR ( ' GREEN ' ) 

TL=. 25 

IF( ICOLOR . EQ . 0) CALL  MRSCOD ( TL , NMRK , RAT) 

IF( Y3 ( 1 ) . NE . 1 . E 100 . AND . ITYPE .EQ.O) CALL  CURVF( X , Y3 , NPT , IMARK ) 
CALL  RESET (' DASH' ) 

IF(  ICOLOR. EQ.DCALL  SETCLR(  '  GREEN '  ) 

IMARK=-1 

IF( N4 . GT . 0 ) CALL  CURVE ( X4 , Y4 , N4 , IMARK) 

CALL  ENDPL(JJ) 

RETURN 

1  FORMAT(6E10.3) 

2  FORMATC  PLOT  SCALES  -  '  , 6E13 . 5) 

END 
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SUBROUTINE  PREDCT(NPLOTP) 

C 

C  THIS  SUBROUTINE  DOES  A  STRAIGHT  LINE  EXTRAPOLATION  OF  DATA  FROM 

SENSORS  2  AND  3  TO  THE  POSITION  OF  SENSOR  1  AND  THEN  THE  DATA 
FROM  THE  PREDICTION  IS  COMPARED  WITH  THE  DATA  OF  SENSOR  1.  NEXT, 

C  THE  MAGNITUDES  AND  PHASE  ANGLES  OF  BOTH  SETS  OF  DATA  ARE  CALCULATED. 

C  IN  THE  VICINITY  OF  LKEY  (A  TIME  OF  INTEREST,  USUALLY  SHOT  EJECTION) 

C  THE  DATA  ARE  AVERAGED  OVER  SHORT  TIME  SPANS  AND  THEIR  MAGNITUDES 

C  AND  PHASE  ANGLES  CALCULATED  AND  PRINTED.  PLOTS  ARE  GENERATED  IF 

C  PLOT?  .NE.  0. 

C 

DIMENSION  AVG(4, 10) , AMPL( 2 , 10) , ANGL( 2 , 10) ,MDAT( 9 ) , NDAT( 9 ) 

COMMON/LEVI /DAT ( 6,2,2500) , LABL( 3,2,24) ,MIN(9) ,NIN(9) ,TSTART, DTIME, 

*  NPT, NPLT, LKEY 

PI2=6. 283 18530 7 17986 5 
RF.AD(5, 1)DZ31  ,DZ32 
WRITE (6 , 2)DZ31 , DZ32 
RAT=DZ31/DZ32 
DO  200  L=1 , NPT 
DO  100  N-1,2 

DAT(4,N,L)=DAT(3,N,L)+RAT*(DAT(2,N,L)-DAT(3,N,L)) 

100  CONTINUE 

DAT (6 , 1 ,L)=SQRT(DAT(1 , l ,L)**2+DAT( 1 , 2 ,L)**2) 

DAT( 5 , l ,L)=S0RT(DAT(4, 1 ,L)**2+DAT(4, 2 ,L)**2) 

IF ( DAT ( 1 , 1 ,L) .EO. .0. AND. DAT (1 , 2 ,L) .EQ. .0)DAT( 1 ,2,0*1 .E-5 
DAT( 6 , 2 ,L)=ATAN2 (DAT( 1 , 1 ,L) , DAT( l , 2 ,L) ) 
IF(DAT(4,l,L).E0..0.AND.DAT(4,2,L).EQ..0)DAT(4,2,L)=l.E-5 
DAT(5 , 2 ,L)=ATAN2(DAT(4, 1 ,L) ,DAT(4, 2 ,L) ) 

IF(DAT(6 , 2 ,L) .LT. .0)DAT(6 , 2 ,L)=DAT(6 , 2 ,L)+PI2 
IF(DAT(5,2,L).LT..0)DAT(5,2,L)=DAT(5,2,L)+PI2 
DAT(6,2,L)=DAT(6,2,L)*360./PI2 
DAT(5,2,L)=DAT(5,2,L)*360./PI2 
200  CONTINUE 
INC=7 

IF( DTIME .EQ .. 02 ) INC=5 

INTV=4 

LINC=INC 

IF( DTIME. EQ.. 02 )LINC=2 
LST=LKEY-INC-( INTV-1 )*LINC+1 
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LSP=LKEY+INTV*LINC 
WRITE( 6 ,4) 

J=0 

DO  400  L=LST,LSP,LINC 

J=J+1 

S1  =  .0 

S2=.0 

S3=.0 

S4=.0 

DO  300  1=1, INC 
S1=S1+DAT( 1 , 1 ,L+I-1 ) 

S2=S2+DAT( 1 , 2 ,L+I-1 ) 

S3=S3+DAT(4, l.L+I-l) 

S4=S4+DAT(4,2,L+I-1) 

300  CONTINUE 

AVG( 1 , j)=Sl/FLOAT( INC) 

AVG( 2 , J)=S2/FL0AT( INC) 

AVG( 3 , J)=S3/FL0AT( INC) 

AVG(4,J)=S4/FL0AT(INC) 

AMPL(1 , J)=SQRT(AVG( 1 , J)**2+AVG(2 , J)**2) 

AMPL(2, J)=SQRT(AVG(3, J)**2+AVG(4, J)**2) 
ANGL(1,J)=ATAN2(AVG(1,J),AVG(2,J)) 

ANGL( 2 , J)=ATAN2 ( AVG( 3 , J) , AVG(4 , J ) ) 

IF(ANGL(1,J).LT..0)ANGL(1,J)=ANGL(1,J)+PI2 

IF(ANGL(2,J).LT..0)ANGL(2,J)=ANGL(2,J)+PI2 

ANGL(1,J)=ANGL(1, J)*360./PI2 

ANGL(2,J)=ANGL(2,J)*360./PI2 

LMID=L+INC/2 

WRITE(6,3)LMID,AVG(1,J),AVG(2,J),AMPL(1,J),ANGL(1,J),AVG(3,J), 
*  AVG(4,J) ,AMPL(2,J) ,ANGL(2,J) 

400  CONTINUE 

IF ( NPLOTP . EQ • 0 ) GO  TO  600 
WRITE(6,7) 

IUPLT=3  ♦ 

READ(IUPLT,6)TIMEL,TIMER 

RE AD ( I UPLT , 5 ) LAX I S , ICOLOR 

NCURV=2 

N4=l 
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DO  500  1=1 ,NPLOTP 

ISC=1 

ITYPE=0 

IF(r.GE.5)lTYPE=l 
IF(I.GE.6)NCURV=1 
NCRV=NCURV*( 1TYPE+1 ) 

RE AD ( IUPLT ,  5  )  ! MDAT ( K ) , NDAT (  K  ) ,K=1 ,NCRV) 

WRITE( 6 ,8) (MDAT(K) , NDAT(K) ,K=1,NCRV) 

NCVV=NCURV 

IF(I.EQ.7)N4=2 

CALL  DATA4 ( NCVV , MDAT , NDAT ,  ISC , ITYPE , LAXIS , ICOLOR, N4 , TIMEL , TIMER , 

*  IUPLT) 

500  CONTINUE 
600  RETURN 

1  FORMAT! 8F1 0.2) 

2  FORMAT!//,'  DZ31 ,  DZ32  \2F10.2) 

3  FORMAT! I5,9X,3F9.5,F12.2,10X,3F9.5,F12.2) 

4  FORMAT!//'  L  BRL :  HORIZONTAL  VERTICAL  MAGNITUDE  ORIENTATION', 

*  '  HEL:  HORIZONTAL  VERTICAL  MAGNITUDE  ORIENTATION') 

5  FORMAT! 2013) 

6  FORMAT! 2F5.0) 

7  FORMAT!//,'  PLOTTING  IN  SUBROUTINE  PPEDCT') 

8  FORMAT!/,'  "MDAT,  NDAT"  PAIRS  -',2013) 

END 
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SUBROUTINE  UNITV(NPLOTU , IPR INT, TST , LST , LSP , IU) 

C 

C  THIS  SUBROUTINE  CALCULATES  THE  TRANSLATION  AND  ROTATION 

C  PARAMETERS  AS  DESCRIBED  IN  SECTION  3.2.3,  EQ.(3.16)  IN 

C  PARTICULAR 

C 

DIMENSION  MDAT (  9  )  , NDAT ( 9 ) 

COMMON/LEVI /DAT( 6, 2, 2500) ,LABL( 3 , 2 , 24) ,MIN(9) , NIN(9 ) ,T START, DTIME , 
*  NPT , NPLT , LKEY 
DO  400  J=1 , 3 
DO  100  L=1 , NPT 

DEN=SQRT(1.+DAT(3,1,L)**2+DAT(3,2,L)**2) 

DN=SQRT(1.+DAT(3,1,L)**2) 

IF( J.GT. 1 )GO  TO  60 

DAT(6,1,L)=1./DN 

DAT(6,2,L)=.0 

DAT( 1 , 1 , L)=-DAT(3 , 1 , L)*DAT( 6 , 1 ,L) 

GO  TO  100 

60  IF ( J.GT. 2) GO  TO  80 
DAT(6,2,L)=DN/DEN 
DAT( 1 , 1 ,L)=-DAT(3,2,L)/(DN*DEN) 

DAT ( 6 , 1 , L ) =DAT ( 1 , l , L ) *DAT ( 3 , 1 , L ) 

GO  TO  100 

80  DAT( 1 , 1 ,L)=1 . /DEN 

DAT ( 6 , 1 , L ) =DAT ( 3 , 1 , L ) *DAT ( 1 , 1 , L ) 

DAT( 6 , 2 ,L)=DAT( 3, 2 ,L)*DAT( 1 , 1 ,L) 

100  CONTINUE 

IF ( NPLOTU .EO.0.OR.J.NE.3) GO  TO  300 
WRITE( 6 , 3) 

IUPLT=10 
REWIND  IUPLT 

READ( iuplt, 2)timel, timer 
READ ( IUPLT , 1 ) LAXI S , ICOLOR 
ISC=0 
NCURV=1 

DO  200  1=1, NPLOTU 
ITYPE=0 

IF( I ,EQ. 3) ITYPE=1 
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LAXIS=0 

IF(I.LE.2)LAXIS=1 

N4=l 

IF( ITYPE .EQ. 1 )N4=2 

NCRV=NCURV*( ITYPE+1 ) 

readCiuplt, i)(mdat(k) ,NDAT(K) ,K=1 ,NCRV) 

WRITE (6 ,4) (MDAT(K) , NDAT(K) ,  K=1 , NCRV) 

ISC=0 

IF( ITYPE.EQ. 1 ) ISC=1 
IF(l.EQ.3.0R.I.EQ.4)lSC=l 

CALL  DATA4(NCURV, MDAT, NDAT, ISC, ITYPE, LAXIS,  ICOLOR,N4,TIMEL, TIMER, 
*  iuplt) 

200  CONTINUE 

300  IF( I PRINT .EQ.O)GO  TO  400 
NARR=3 
MDAT( 1  )=6 
MDAT(2)=6 
MDAT ( 3 ) = 1 
NDAT ( 1 ) = 1 
NDAT(2)=2 
NDAT(3)=1 

CALL  WRFILE(OAT,MDAT,NDAT,NARR,TST,LST,LSP,LABL,OTIME,IU) 

400  CONTINUE 
500  RETURN 

1  FORMAT (2013) 

2  FORMAT( 2F5 . 0 ) 

3  FORMAT(//,'  PLOTTING  IN  SUBROUTINE  UNITV') 

4  FORMAT(/,'  "MDAT,  NDAT"  PAIRS  -',2013) 

END 
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SUBROUTINE  WRFILE ( DAT ,MDAT , NDAT , NARR , TST , LST , LSP , LABL , DTIME , III) 

C 

C  THIS  SUBROUTINE  WRITES  DATA  ARRAYS  TO  AN  OUTPUT  FILE  TO  BE  SAVED 

C  FOR  FUTURE  USE 

C 

C  NARR  =  NUMBER  OF  ARRAYS  TO  BE  WRITTEN 

C  MDAT  AND  NDAT  CONTROL  THE  ORDER  IN  WHICH  THE  ARRAY  IS  WRITTEN 

C  MDAT  =  FIRST  INDEX  OF  3-DIMENSIONAL  ARRAY  DAT 

C  NDAT  =  SECOND  INDEX  OF  3-DIMENSIONAL  ARRAY  DAT 

C  LST  AND  LSP  CONTROL  THE  AMOUNT  OF  DATA  WRITTEN  BY  SPECIFING  THE 

C  STARTING  AND  STOPPING  ELEMENTS  OF  THE  THIRD  INDEX  OF  THE 

C  3-DIMENSIONAL  ARRAY  DAT 

C  LST  *  STARTING  ELEMENT 

C  LSP  =  STOPPING  ELEMENT 

C 

DIMENSION  DAT( 6, 2,2500) ,MDAT(9) ,NDAT(9) ,LABL(3 ,2,24) ,X(500) 

DO  200  1=1, NARR 
K=0 

DO  100  L=LST, LSP 
K=K+1 

X(K)=DAT(MDAT(I),NDAT(I) ,L) 

100  CONTINUE 

WRITE ( 10 )( LABL ( 1 , 1 , J) , J=1 ,24) 

write( iu)k 

WRITE ( IU) TST, DTIME 
WRITE(IU)(X(J),J=1,K) 

WRITE ( 6,1 )MDAT( I) , NDAT(l) , (LABL( 1 , 1 , j) , J=1 , 24) 

200  CONTINUE 
RETURN 

1  FORMAT ('  WRITING  ARRAY  ON  FILE  -  INDICES  ARE  ' , 213 , 5X, 24A2) 

END 


C.2.3.  Listing  Of  Auxiliary  FORTRAN  Subroutines  (MFA  File  DIFFATR) 


1 

C 

C 

c 

5  C 

C 
C 


10 


15 


co 


20 


SUBROUTINE  DIFFER (Z , N , DTIME , IUD , IUF ) 

IF  IUD  .NE.  0,  THIS  SUBROUTINE  DIFFERENTIATES  AN  ARRAY  USING  AN 
FIR  DIFFERENTIATOR;  USUALLY,  THE  ONE  IN  FILE  DIFFl  IS  USED. 


IF  IUF  .NE.  0,  THEN  THE  ARRAY  IS  LOWPASS  FILTERED. 


DIMENSION  Z(N),H(512) 

PI2=6. 283185307179865 
IF ( IUD . EQ . 0 )GO  TO  150 
C=PI2 /DTIME 

CALL  COEF2(NFLT,H, IUD, JTPE) 
REWIND  IUD 
DO  100  1=1 ,NFLT 
H(I)=H(I)*C 
100  CONTINUE 

CALL  CONVOL(H, Z ,N, NFLT , JTPE ) 
150  IF ( IUF .EQ . 0)GO  TO  200 
CALL  FILT(Z ,N, IUF ) 

200  RETURN 
END 
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SUBROUTINE  C0EF1 (NFILT, H , IU) 

READS  IN  COEFFICIENTS  OF  AN  FIR  LOWPASS  OR  HIGHPASS  FILTER 

DIMENSION  H(512) 

READ ( IU) NFILT , FP , TBWID , DP , DS 
IF( ISW. NE . 1 ) WRITE ( 6,1) 

IF(ISW.NE. 1)WRITE(6, 2)NFILT,FP, TBWID, DP, DS 
N=(NFILT+l)/2 

read(iu)(h(i),i=i,n) 

IF( ISW.NE. 1)WRITE(6 , 3) 

IF(ISW.EQ.1)G0  TO  200 
DO  100  1=1, N, 5 
IST=I+4 

IF(IST.GT.N)IST=N 

WRITE ( 6, 4) I, (H( II) , I 1=1, 1ST) 

100  CONTINUE 
WRITE(6,5) 

200  ISW=1 
RETURN 

1  FORMAT (// , '  FILTER  PARAMETERS,') 

2  FORMAT( 5X, '  NFILT  =  ',15,'  FP  =  '.F10.5,'  TBWID  =  '.F10.5, 
*'  DP  =  ' ,F10.5, '  DS  =  ' ,F10.5) 

3  FORMAT( 5X, '  FILTER  COEFFICIENTS') 

4  FORMAT (5X,I5,5E15.8) 

5  FORMAT(//) 

END 


1 

5 

10 

15 

20 

25 

30 

35 


O  CJ 


SUBROUTINE  C0EF2 ( NFILT , H , IU , JTYPE ) 


C 


READ  IN  COEFFICIENTS  OF  AN  FIR  DIFFERENTIATOR 


DIMENSION  H(64) , EDGE (20) ,FX( 10) ,WTX( 10) 

IF ( ISW.NE . 1 ) WRITE ( 6,1) 

READ( IU) NFILT , JTYPE , NBANDS 

I F ( I SW . NE . 1 ) WRITE ( 6 , 2 ) NF I LT , JTYPE , NB AND  S 

JB=2*NBANDS 

READ(IU)(EDGF,(J)  ,J=1 ,JB) 

IF ( ISW.NE . 1 ) WRITE ( 6 , 5 ) (EDGE( J) ,J=1,JB) 

READ(IU)(FX(J) ,J=1, NBANDS) 

IF( ISW.NE . 1 ) WRITE ( 6 ,6)(FX(j) ,J=1, NBANDS) 

RE AD ( IU ) ( WTX ( J ) , J= 1 , NBANDS ) 

I F ( I SW . NE . 1 ) WRITE ( 6 , 7 ) ( WTX ( J ) , J= 1 , NBANDS ) 

READ(lU)(H(j) , J=1 ,64) 

IF ( ISW.NE . 1 ) WRITE (6 ,3) 

IF( ISW.EQ. l)GO  TO  200 
N=NFILT/ 2+1 
DO  100  1*1, N, 5 
IST=I+4 

IF(IST.GT.N)IST=N 
WRITE(6,4)I,(H(II) ,II=I,IST) 

100  CONTINUE 
WRITE (6, 8) 

200  ISW=1 
RETURN 

1  FORMAT(//,'  FILTER  PARAMETERS') 

2  FORMAT( 5X, '  NFILT  =  ',15,'  JTYPE  =  ',15,'  NBANDS  =  ’,15) 

3  F0RMAT(5X, '  FILTER  COEFFICIENTS') 

4  FORMAT ( 5X,I5,5E15.8) 

5  FORMAT(5X,'  EDGE  '.20F6.4) 

6  FORMATC 5X, '  FX  ’,10F6.2) 

7  FORMAT (5X, '  WTX  * ,10F7.1) 

8  FORMATC//) 

END 


1 

5 

10 

15 

20 

25 

30 

35 

40 


SUBROUTINE  CONVOL( H , X , NDIM , NFILT , JTYPE ) 

C 

C  THIS  CONVOLUTION  IS  VALID  ONLY  FOR  NFILT  =  ODD  INTEGER. 

C  NFILT  MAX.  SET  TO  1023  -  CAN  BE  RESET  BE  REDIMENBIONING  ARRAYS 

C  IF  JTYPE  =  1,  EVEN  SYMMETRY  ASSUMED  FOP  REFLECTION  AT  ENDPOINTS 

C  IF  JTYPE  =  2,  ODD  SYMMETRY  ASSUMED  FOR  REFLECTION  AT  ENDPOINTS 

C 

DIMENSION  H(512),X(NDIM) ,S( 1030), T( 1030) 

IF( JTYPE. EQ. 1 )SGN=1 . 

IF(JTYPE.EQ.2)SGN=-1 . 

IF ( JTYPE .EQ. 1 )A=0 . 

IF( JTYPE.EQ. 2)A=2 . 

IF( JTYPE. LT.l .OR. JTYPE. GT. 2)WRITE(6, 1) 

IF(JTYPE.LT. 1 .OR. JTYPE. GT. 2) STOP 

L=NFILT-1 

NCOEF=NFILT/2+l 

J=NCOEF 

K=J-1 

DO  10  1=1 ,NCOEF 
S(l)=SGN*X(NCOEF-I+l)+A*X(l) 

S(NFILT-I+1)=X(J) 

T(NCOEF+I-l)=SGN*X(NDIM-I+l)+A*X(NDIM) 

J=J-1 

10  CONTINUE 

DO  40  1=1, NDIM 

X(I)=.0 

DO  20  J=1,K 

X( I  )=X( l)+H( J)*( SGN*S( J)+S (NFILT— J+l ) ) 

20  CONTINUE 

X(l)=X(l)+H(NCOEF)*S(NCOEF) 

IF ( I. EO. NDIM) GO  TO  50 
DO  30  J=1 , L 
S(J)=S(J+1) 

30  CONTINUE 

IF(l.LE . NDIM-NCOEF) s (NFILT)=X( i+ncoef) 

IF( I . GT . NDIM-NCOEF) S(NFILT)=T( I+NFILT-NDIM+1 ) 

40  CONTINUE 
50  RETURN 

1  FORMAT ('  ERROR  IN  SUB  CONVOL  -  JTYPE  NOT  EQUAL  TO  1  OR  2') 

END 


n  o 


1  SUBROUTINE  FILT(X,NPT, IU)  ’ 

THIS  SUBROUTINE  LOWPASS  OR  HIGHPASS  FILTERS  A  DATA  ARRAY 
C 

5  DIMENSION  H(512) ,X(NPT) 

JTYP=1 

CALL  COEF1 ( NFILT,H, IU) 

REWIND  IU 

200  CALL  CONVOL ( H , X , NPT , NF ILT , JTYP ) 

10  RETURN 

END 


CO 

PO 


C.2.4.  Card  Image  Formats  That  Are  Required  For  The  Input  Data  Section  of  MFA  File  MUZPRED/UN= BOOTS 

TABLE  C.2.  INPUT  DATA  FORMATS  FOR  TAPE  UNIT  5 


Card 

Condition* 

Column 

Format 

Variable 

Calling  Subroutine 

Description 

i 

1-3 

13 

NCH 

MAIN 

No.  of  data  sets  to  be  inputted 

4-6 

13 

NPLOTP 

No.  of  plots  to  be  generated  in 
sub.  PREDCT 

7-9 

13 

NPLOTC 

No.  of  plots  to  be  generated  in 
sub.  CALCZ 

10-12 

13 

NPLOTU 

No.  of  plots  to  be  generated  in 
sub.  UNITV 

13-15 

13 

IPRINT 

If  IPRINT  f09  output  data  file 
is  generated. 

2 

1-3 

13 

MIN(I) 

MAIN 

First  index  for  placement  of 
data  into  array  DAT 

4-6 

13 

NIN(I) 

Second  index  for  placement  of 
data  into  array  DAT 

7-9 

13 

IU 

Tape  unit  no.  on  which  data 
file  resides 

3 

Repeat  card 

2  until  NCH 

cards  are 

inputted . 

4 

1-5 

15 

LKEY 

MAIN 

Index  at  which  to  make  time  =  0. 

5 

1-10 

F10 .2 

DZ31 

PREDCT 

Distance  from  sensor  3  to  sensor 
1 ,  cm 

11-20 

CM 

O 

DZ32 

PREDCT 

Distance  from  sensor  3  to  sensor 
2 ,  cm 

6 

1-5 

15 

LZM 

CALCZ 

Index  of  start  of  recoil  motion 
(earlier  experiment) 

6-10 

15 

LPM 

Index  of  peak  pressure 
(earlier  experiment) 

11-15 

15 

LEM 

Index  of  shot  exit 

(earlier  experiment) 

*  Omit  card  if  condition  is  not  met. 


TABLE  C.2.  INPUT  DATA  FORMATS  FOR  TAPE  UNIT  5  (continued) 


co 

-p* 


Card  Condition 

Column 

Format 

Variable 

Calling  Subroutine 

Description 

7 

1-5 

15 

LPD 

CALCZ 

Index  of  peak  pressure 
(this  experiment) 

6-10 

15 

LED 

Index  of  shot  exit 

- 

(this  experiment) 

8 

1-5 

15 

MPOLY 

CALCZ 

Degree  of  polynomial  to  be 
inputted  (earlier  experiment) 

9 

1-10 

E10.3 

C(l) 

CALCZ 

Coefficients  of  polynomial 
(earlier  experiment) 

11-20 

E10.3 

C  (  2) 

• 

• 

• 

• 

• 

• 

• 

• 

s 

• 

0 

C(MP0LY+1) 

10 

1-10 

E10.3 

Z (1 , 1 ) 

CALCZ 

Distance  of  sensor  1  (horiz)  fro 
muzzle  face,  cm 

11-20 

E10.3 

Z (1 , 2) 

Distance  of  sensor  l(vert.) 
from  muzzle  face,  cm 

21-30 

E10.3 

Z(2, 1) 

Distance  of  sensor  2  from 
muzzle  face,  cm 

31-40 

E10.3 

Z(3, 1) 

Distance  of  sensor  3  from 
muzzle  face,  cm 

11 

1-3 

13 

NZREF 

CALCZ 

No.  of  times  the  displacement 
and  muzzle-pointing  calculations 
are  to  be  done 

12 

1-10 

E10.3 

ZREF 

CALCZ 

Distance  from  muzzle  face  at 

which  the  displacement  and 
muzzle-pointing  calculations 
are  to  be  done 


TABLE  C.2.  INPUT  DATA  FORMATS  FOR  TAPE  UNIT  5  (continued) 


Card 

Condition 

Column 

Format 

Variable 

Calling  Subroutine 

Description 

13 

IPRINT^O 

1-3 

13 

MDAT(l) 

CALCZ 

First  index  in  array  DAT 
of  first  set  of  data  to  be 

written  to  output  file 

4-6 

13 

NDAT(l) 

Second  index  in  array  DAT 
of  first  set  of  data  to  be 
written  to  output  file 

7-9 

13 

MDAT(2) 

First  index  of  second  set 

10-12 

13 

NDAT ( 2 ) 

Second  index  of  second  set 

13-15 

13 

MDAT ( 3 ) 

First  index  of  third  set 

16-18 

13 

NDAT (3) 

Second  index  of  third  set 

14  Repeat  cards  12  and  13  for  NZREF  times. 


oo 

cn 


TABLE  C.3.  INPUT  DATA  FORMATS  FOR  TAPE  UNITS  3,7,  AND  10 


Card 

Condition* 

Column 

Format 

Variable 

Description 

i 

1-5 

F5.0 

TIMEL 

Time  to  start  plotting  (time  =  .0  at  LKEY),  ms 

6-10 

F5.0 

TIMER 

Time  to  stop  plotting,  ms 

2 

1-3 

13 

LAXIS 

=1  Plot  axis  labels  to  be  inputted 

Plot  axis  labels  are  blank 

4-6 

13 

ICOLOR 

=0  Black/white  plot 
=1  Color  plot 

3 

1-3 

13 

MDAT(l) 

First  index  in  array  DAT  of  first  set  of  data  to 

be  plotted 

4-6 

13 

NDAT ( 1 ) 

Second  index  in  array  DAT  of  first  set  of  data  to 

be  plotted 

7-9 

13 

MDAT ( 2 ) 

First  index  of  second  curve 

10-12 

13 

NDAT (2) 

Second  index  of  second  curve 

s 


00  9 

MDAT(NCRV) 

NDAT(NCRV) 


4 

N4>0 

1-10 

F10.0 

X4(  1 ) 

X-coordinate  of 

a  point  to  be  symbol-plotted 

11-20 

F10.0 

Y4(l) 

Y-coordinate  of 

a  point  to  be  symbol-plotted 

5 

LAXIS=1 

1-10 

A10 

LXNAME 

Label  for  x-axis 

11-20 

A10 

LYNAME 

Label  for  y-axis 

6 

ISC#) 

1-10 

E10.3 

XORIG 

Minimum  value  of 

x-axis 

11-20 

E10.3 

XMAX 

Maximum  value  of 

x-axis 

21-30 

E10.3 

XSTP 

Interval  between 

tic-marks  on  x-axis 

31-40 

E10.3 

YORIG 

Minimum  value  of 

y-axis 

41-50 

E10.3 

YMAX 

Maximum  value  of 

y-axis 

51-60 

E10.3 

YSTP 

Interval  between 

tic-marks  on  y-axis 

7 

Repeat 

cards  3 

through  6 

for  NPLOTP,  NPLOTC,  . 

or  NPLOTU  times,  depending 

on  whether  the  calling  subroutine  is  PREDCT,  CALCZ,  OR  UNITV. 


*  Omit  card  if  condition  is  not  met. 
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APPENDIX  D 


SAMPLE  PROBLEM 

As  alluded  to  earlier  in  this  report,  the  computer  programs  have  been  developed  for  the  analysis  of 
muzzle  displacements  collected  in  a  collaborative  experiment  between  BFL  and  RARDE  in  May  and  June  1984  at 
RARDE,  Fort  Halstead  [5].  This  experiment,  carried  out  with  a  40mm  rifled  gun,  had  the  objective  to  assess 
the  SDT  measurement  technique  as  a  means  to  define  projectile  launch.  It  included  the  measurement  of  the 
projectile  dynamics  with  respect  to  the  muzzle  during  shot  exit  as  well  as  the  motion  of  the  muzzle  with 
respect  to  the  ground.  The  latter  was  monitored  by  three  orthogonal  displacement  transducers  employing  a 
SDT  instrumentation  for  sensor  1  and  electrooptical  displacement  transducers  [A-8]  for  sensors  2  and  3.  The 
three  data  sets  for  round  #5  plus  the  muzzle  recoil  velocity  recorded  by  an  earlier  experiment  (see  Section 
4.3.1)  form  the  input  to  the  sample  problem.  Since  the  result  of  the  analysis  of  the  data  from  the 
collaborative  experiment  [6]  will  be  discussed  in  a  separate  BRL  technical  report,  only  sample  outputs 
including  plots  are  shown  for  this  round. 


oo 

CO 


D.l.  Program  MDZM040 


D.1.1.  Sample  Input  Card  Images 


co 

o 


0  0  0  1 
6  10  1 
1  1  1 
1  2  2 
2  1  8 

2  2  9 

3  1  10 

3  2  11 


202.52  .00221  .00444  215.11 

108.9  242.1 

263.9  227.8 

1  1  .0024414  .9985  1. 

1  2  .0024414  .99925  1. 

2  1  .0048828  .99925  1. 

2  2  .0048828  .99925  1. 

3  1  .0048828  .9985  1. 

3  2  .0048828  1.001  1. 


.00349 


00221 


APPENDIX  D 


D.1.2.  Sample  Printed  Output  (Program  MUZM040) 


IPLT1  =  0  IPLT2  =  0  IPLT3  =  0  IPRT  =  1 

NCH  =  6  ICVT  =  1  IFILT  =  0  IDERIV  =  1 


INPUT  DATA  CHANNELS 

1  1  1 
1  1  BR0512 

4095  0.00000 

1  2  2 

1  2  BR0511 

4095  0.00000 

2  1  8 

2  1  BR0509 

4095  0.00000 

°  2  2  9 

2  2  BR.0508 

4095  0.00000 

3  1  10 

3  1  BR0510 

4095  0.00000 

3  2  11 

3  2  BR0507 

4095  0.00000 


1  1  1BRL  TUBE  HOR 

.01000 

1  1  1BRL  TUBE  VERT 

.01000 

1  1  1HEL  TUBE  FWD  HOP 

.01000 

1  1  1HEL  TUBE  FWD  VERT 

.01000 

1  1  1HEL  TUBE  REAR  HOP 

.01000 

1  1  1HEL  TUBE  REAR  VERT 

.01000 


CALIBRATION  CONSTANTS 

202.52000  .00221  .00444  215.11000  .00349  .00221 

108.90000  242.10000 

263.90000  227.80000 


1 

1 

.00244 

.99850 

1.00000 

.00245 

1 

2 

.00244 

.99925 

1.00000 

.00244 

2 

1 

.00488 

.99925 

1.00000 

.00489 

2 

2 

.00488 

.99925 

1.00000 

.00489 

3 

1 

.00488 

.99850 

1.00000 

.00489 

3 

2 

.00488 

1.00100 

1.00000 

.00488 

1 

1 

.00245 

202.52000 

1 

2 

.00244 

215.11000 

2 

1 

.00489 

108.90000 

2 

2 

.00489 

242.10000 

3 

1 

.00489 

263.90000 

3 

2 

.00488 

227.80000 

FILTER  PARAMETERS 

NFILT  =  31  JTYPE  =  2  NBANDS  =  1 

EDGE  0.0000  .4000 

FX  1.00 

WTX  1 . 0 

FILTER  COEFFICIENTS 

1  . 16166620E-04  -. 74131438E-04  .21102257E-03  -.49246632E-03 

6  -. 18862205E-02  . 32833366E-02  -. 54084061E-02  . 85345464E-02 

11  . 19559239F-01  -. 29218184E-01  .44647104E-01  -. 73733156E-01 

16  0. 


. 10096462E-02 
.  13047 140E-01 
.  15615738E+00 


CHANNEL  BR0512 


1  1  1BRL  TUBE  HOR 


HAS  BEEN  DIFFERENTIATED 


FILTER  PARAMETERS, 


NFILT 

=  127  FP  = 

.00500  TBWID 

FILTER 

COEFFICIENTS 

1  - 

. 5761 1818E-02 

-. 13120665E-02 

6  - 

. 18272702E-02 

-.19407P50E-02  • 

11  - 

. 22685610E-02 

-. 23074309E-02  ■ 

16  - 

. 22201999E-02 

-. 21269670E-02  • 

21  - 

. 14040333E-02 

-. 1 1271255E— 02  ■ 

26 

. 39177683E-03 

.880901 17E-03 

31 

. 32498950E-02 

. 39438586E-02 

36 

. 70536956E-02 

. 79036343E-02 

41 

. 1 1502901E-01 

. 12431251E-01 

46 

. 16127318E-01 

. 17023636E-01 

51 

. 20359358E-01 

. 21 108628E-01 

56 

. 23637514E-01 

. 24135636E-01 

61 

. 25501 145E-01 

. 25673888E-01 

.01541  DP  = 


14458869E-02 

20444951E-02 

23235239E-02 

20019212E-02 

80807552E-03 

14120564E-02 

46685324E-02 

87797060E-02 

13361246E-01 

17900295E-01 

21814687E-01 

24573834E-01 

25780396E-01 


.01000 

-. 1 5779811E-02 
-. 21345667E-02 
-. 23165943E-02 
-. 18400413E-02 
-.45054632E-03 
. 19883328E-02 
. 54340488E-02 
. 96719531E-02 
. 14291479E-01 
.  18750877E-01 
.  22474542E-01 
. 24947050E-01 
.25814786E-01 


DS  =  .01000 

-. 17054889E-02 
- . 221 10148E-02 
-. 22818637E-02 
16423613E-02 
-.49316320E-04 
.  25922726E-02. 

. 62277021E-02 
.  10582451E-01 
.  1 5214561E-01 
.  19572419E-01 
.  23082659E-01 
.  25256526E-01 


CHANNEL 

BR0512 

1 

1 

1BRL 

TUBE 

CHANNEL 

8R0512 

1 

1 

1BRL 

TUBE 

CHANNEL 

BROS  12 

1 

1 

1 RRL 

TUBE 

CHANNEL 

BR0511 

1 

1 

1BPL 

TUBE 

CHANNEL 

BR0511 

1 

1 

1RRL 

TUBE 

CHANNEL 

BR0511 

1 

1 

1BPL 

TUBE 

CHANNEL 

BR0511 

1 

1 

1BRL 

TUBE 

CHANNEL 

BR0509 

1 

1 

1HEL 

TUBE 

CHANNEL 

BR0509 

1 

1 

1HEL 

TUBE 

CHANNEL 

BR0509 

1 

1 

1HEL 

TUBE 

CHANNEL 

BR0509 

1 

1 

1HFL 

TUBE 

CHANNEL 

BR0508 

1 

1 

1HEL 

TUBE 

CHANNEL 

BR0508 

1 

1 

1 HEL 

TUBE 

CHANNEL 

BR0508 

1 

1 

1FEL 

TUBE 

CHANNEL 

BR0508 

1 

1 

1HFL 

TUBE 

CHANNEL 

BR0510 

1 

1 

1HEL 

TUBE 

CHANNEL 

BR0510 

1 

1 

1HEL 

TUBE 

HOR 

HAS 

BEEN 

FILTERED 

HOR 

FAS 

BEEN 

DIFFERENTIATED 

HOP 

HAS 

BEEN 

FILTERED 

VERT 

HAS 

BEEN 

DIFFERENTIATED 

VFRT 

HAS 

BEEN 

FILTFRED 

VERT 

HAS 

BEEN 

DIFFERENTIATED 

VERT 

HAS 

BEEN 

FILTERED 

FWD  HOR 

HAS 

BEEN 

DIFFERENTIATED 

FWD  HOR 

HAS 

BEFN 

FILTERED 

FWD  HOP 

HAS 

BEEN 

differentiated 

FWD  HOR 

HAS 

BEEN 

FILTERED 

FWD  VERT 

HAS 

BEEN 

DIFFERENTIATED 

FWD  VERT 

HAS 

BEEN 

FILTERED 

FWD  VERT 

HAS 

BEEN 

DIFFERENTIATED 

FWD  VERT 

HAS 

BEEN 

FILTERED 

REAR  HOR 

HAS 

BEEN 

DIFFERENTIATED 

REAR  HOR 

HAS 

BEEN 

FILTERED 

o 

CO 


CHANNEL 

BR0510 

1 

1 

1HEL 

TUBE 

REAR 

HOR 

HAS 

BEEN 

CHANNEL 

BR0510 

1 

1 

1HEL 

TUBE 

REAR 

HOR 

HAS 

BEEN 

CHANNEL 

BR0507 

1 

1 

1HEL 

TUBE 

REAR 

VERT 

HAS 

BEEN 

CHANNEL 

BR0507 

1 

1 

1HEL 

TUBE 

REAR 

VERT 

HAS 

BEEN 

CHANNEL 

BP0507 

1 

1 

1HEL 

TUBE 

REAR 

VERT 

HAS 

BEEN 

CHANNEL 

BR0507 

1 

1 

1HEL 

TUBE 

REAR 

VERT 

HAS 

BEEN 

INPUT  DATA  -  UNFILTERED 


WRITING 

ON 

FILE 

- 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

WRITING 

ON 

FILE 

- 

BR0511 

1 

1 

1 BRL 

TUBE 

VERT 

WRITING 

ON 

FILE 

- 

BR0509 

1 

1 

1HEL 

TUBF 

FWD  FOR 

WRITING 

ON 

FILE 

- 

BR0508 

1 

1 

1 HFL 

TUBE 

FWD  VERT 

WRITING 

ON 

FILE 

- 

BR0510 

1 

1 

1HEL 

TUBE 

REAR  HOR 

WRITING  ON  FILE  -  BR0507 

FIRST  DERIVATIVE  -  FILTERED 

l 

1 

1HEL 

TUBE 

REAR  VERT 

WRITING 

ON 

FILE 

- 

BR0512 

1 

1 

1BPL 

TUBE 

HOR 

WRITING 

ON 

FILE 

- 

BR0511 

1 

1 

1  BRL 

TUBE 

VERT 

WRITING 

ON 

FILE 

- 

BR0509 

1 

1 

1HEL 

TUBE 

FWD  HOF 

WRITING 

ON 

FILE 

- 

BR0508 

1 

1 

1HEL 

TUBE 

FWD  VERT 

WRITING 

ON 

FILE 

- 

BR0510 

1 

1 

1HEL 

TUBE 

REAR  HOR 

WRITING  ON  FILE  -  BR0507 

SECOND  DERIVATIVE  -  FILTERED 

1 

1 

1HEL 

TUBE 

REAR  VERT 

WRITING 

ON 

FILE 

- 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

WRITING 

ON 

FILE 

- 

BR0511 

1 

1 

1BRL 

TUBE 

VERT 

WRITING 

ON 

FILF 

- 

BR0509 

1 

1 

1HEL 

TUBE 

FWD  HOR 

WRITING 

ON 

FILE 

- 

BR0508 

1 

1 

1HEL 

TUBE 

FWD  VERT 

WRITING 

ON 

FILE 

- 

BR0510 

1 

1 

1HEL 

TUBE 

REAR  HOR 

WRITING 

ON 

FILE 

- 

BR0507 

1 

1 

1HEL 

TUBE 

REAR  VERT 

DIFFERENTIATED 

FILTERED 

DIFFERENTIATED 

FILTERED 

DIFFERENTIATED 

FILTERED 


END  OF  DISSPLA  8.2  --  0  VECTORS  GENERATED  IN  0  PLOT  FRAMES. 
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APPENDIX  D 


D.2.  Program  MUZPRED 


D.2.1.  Sample  Input  Card  Images  From  Tape  Unit  5  (Normal  Input) 


10 


8  7  6 

1  1  1 
1  2  1 
2  1  1 
2  2  1 
3  1  1 

3  2  1 

1  1  4 

1  2  4 

2059 

-45.36 
2175  2342 
1490  2059 
5 

33679E-2 

-2.94 

2 

.0 

6  2  2 
-7.52 
6  2  2 


3  1 

-12.00 
3474 

3.7657 
-3.25 

2  2 
2  2 


.005 


.41978 

-36.56 


-15.709 

-48.56 


23.535  -10.359 


D.2.2  Sample  Input  Card  Images  From  Tape  Unit  3  (MFA  File  DISPL40) 


<£> 


-9.  3. 

1  0 

114  1 

2059. 
TIME,  MS 
-9. 

12  4  2 

2059. 
TIME,  MS 
-9. 

6  15  1 

2059. 
TIME,  MS 
-9. 

6  2  5  2 

2059. 
TIME,  MS 
-9. 

1112 

2059. 

HORIZONTAL 

-.08 

1112 

2059. 

HORIZONTAL 

-.08 

1112 

2059. 

HORIZONTAL 

-.02 


.0 

HORIZONTAL 

3.  3. 

.0 

VERTICAL 
3.  3. 

.05 

MAGNITUDE 

3.  3. 

260. 

ANGLE 

3.  3. 

4  14  2 

.0 

VERTICAL 
.08  .04 

.0 

VERTICAL 
.08  .04 

.0 

VERTICAL 
.03  .01 


-.06 


-.10 


.0 


60. 


-.10 


-.10 


-.02 


.04 


.06 


.10 


360. 


.06 


.06 


.02 


.04 


.02 


100. 


.04 


.04 


03 


01 


D.2.3.  Sample  Input  Card  Images  From  Tape  Unit  7  (MFA  File  RECL40) 


-9.  3. 

1  0 
2  1 
2059. 
(0)MSC1) 

2  2 
2059. 
(0)MS(1) 

2  12  2 
2059. 
(0)CM(1) 
-.02 

3  1 

2059. 
(0)MS(1) 

-10. 

3  2 

2059. 
(0)MS(1) 

3  13  2 

2059. 

(0)CM/MS(1) 

-.02 


.0 

(0)CM(1) 

.0 

(0)CM(1) 


.0 

(0)CM(1) 

.02 

.01 

-.02 

.02 

.0 

(0) cm/msC 1) 

4. 

2. 

-.015 

.015 

.0 

(0)CM/MS( 1) 

.0 

(0) CM/MS( 1 ) 

.02 

.01 

-.02 

.02 

.01 


.005 


.01 


D.2.4.  Sample  Input  Card  Images  From  Tape  Unit  10  (MFA  File  UVECT40) 


-1.  1. 

1  0 
6  1 

2059.  .0 

(0)MS(1) 

6  2 

2059.  .0 

(0)MS(1) 

6  16  2 

2059.  .0 

-.0005  .0005 


0005 


-.0005 


0005 


0005 


APPENDIX  D 


D.2.5.  Sample  Printed  Output  (Program  MUZPRED) 

NCH,  NPLOTP,  NPLOTC,  NPLOTU,  IPRINT  -87631 


MIN,  NIN,  IU  - 

1  1 

1 

1  1  BR0512 

i 

i 

1 BRL 

TUBE  HOE 

NDATA ,  TSTAPT, 

DTIME 

- 

4095 

0.00000 

.01000 

MIN,  NIN,  IU  - 

1  2 

1 

1  2  BR0511 

i 

i 

1BRL 

TUBE  VERT 

NDATA,  TSTART, 

DTIME 

- 

4095 

0.00000 

.01000 

MIN,  NIN,  IU  - 

2  1 

1 

2  1  BR0509 

i 

i 

1HEL 

TUBE  FWD 

HOR 

NDATA,  TSTART, 

DTIME 

- 

4095 

0.00000 

.01000 

MIN,  NIN,  IU  - 

2  2 

1 

2  2  BR0508 

i 

i 

1HEL 

TUBE  FWD  1 

VERT 

NDATA,  TSTART, 

DTIME 

- 

4095 

0.00000 

.01000 

MIN,  NIN,  IU  - 

3  1 

1 

3  1  BR0510 

i 

i 

IHEL 

TUBE  REAR 

HOR 

NDATA,  TSTART, 

DTIME 

- 

4095 

0.00000 

.01000 

MIN,  NIN,  IU  - 

3  2 

1 

3  2  BR0507 

i 

i 

IHEL 

TUBE  REAR 

VERT 

NDATA,  TSTART, 

DTIME 

- 

4095 

0.00000 

.01000 

MIN,  NIN,  IU  - 

1  1 

4 

1  1  BR0512 

i 

i 

1BRL 

TUBE  HOR 

NDATA,  TSTART, 

DTIME 

- 

4095 

0.00000 

.01000 

MIN,  NIN,  IU  - 

1  2 

4 

1  2  BR0511 

i 

i 

1BRL 

TUBE  VERT 

NDATA,  TSTART, 

DTIME 

- 

4095 

0.00000 

.01000 

LKEY  =  2059 


-45.36  -12.00 


DZ31 ,  DZ32 


L 

BRL : 

HORIZONTAL 

VERTICAL 

MAGNITUDE 

ORIENTATION 

2035 

-.02361 

-.06044 

.06489 

201.34 

2042 

-.02248 

-.06166 

.06563 

200.03 

2049 

-.02118 

-.06319 

.06664 

198.53 

2056 

-.02036 

-.06520 

.06830 

197.34 

2063 

-.02049 

-.06730 

.07035 

196.93 

2070 

-.02152 

-.06952 

.07278 

197.20 

2077 

-.02317 

-.07164 

.07530 

197.92 

2084 

-.02448 

-.07340 

.07737 

198.44 

PLOTTING  IN 

SUBROUTINE 

PREDCT 

"MDAT, 

NDAT" 

PAIRS  -  1 

1  4  1 

SPECIAL  COORDINATES  20.590  0.000 

PLOT  LIMITS  -  1160  2359  -9.0000  2.9900 

TIME,  MS  HORIZONTAL 

PLOT  SCALES  -  -.90000E+01  .30000E+01  .30000E+01 

"MDAT,  NDAT"  PAIRS  -  1  2  4  2 

SPECIAL  COORDINATES  20.590  0.000 

PLOT  LIMITS  -  1160  2359  -9.0000  2.9900 

TIME,  MS  VERTICAL 

PLOT  SCALES  -  -.90000E+01  .30000E+01  .30000E+01 

"MDAT,  NDAT"  PAIRS  -  6  1  5  1 

SPECIAL  COORDINATES  20.590  .050 

PLOT  LIMITS  -  1160  2359  -9.0000  2.9900 

TIME,  MS  MAGNITUDE 

PLOT  SCALES  -  -.90000E+01  .30000E+01  .30000E+01 

"MDAT,  NDAT"  PAIRS  -  6  2  5  2 

SPECIAL  COORDINATES  20.590  260.000 

PLOT  LIMITS  -  1160  2359  -9.0000  2.9900 

TIME ,  MS  ANGLE 

PLOT  SCALES  -  -.90000E+01  .30000E+01  .30000E+01 


HORIZONTAL 

VERTICAL 

MAGNITUDE 

ORIENTATION 

-.01174 

-.04890 

.05029 

193.51 

-.01166 

-.05164 

.05294 

192.73 

-.01212 

-.05455 

.05588 

192.52 

-.01292 

-.05751 

.05894 

192.66 

-.01385 

-.06032 

.06180 

192.93 

-.01472 

-.06288 

.06458 

193.17 

-.01546 

-.06506 

.06687 

193.37 

-.01616 

-.06666 

.06859 

193.63 

-.60000E-01  . 20000E-01  .40000E-01 


-.  10000E+00  .40000E-01  .60000F-01 


0.  . 20000E-01  . 10000E+00 


60000E+02 


10000E+03 


36000E+03 


"MDAT,  NDAT"  PAIRS  -  1  1  1  2  4  1  4  2 

SPECIAL  COORDINATES  20.590  0.000 

PLOT  LIMITS  -  1160  2359  -.0236  -.0291 

HORIZONTAL  VERTICAL 

PLOT  SCALES  -  -.80000E-01  .40000E-01  .80000E-01  -.10000E+00  .40000E-01  .60000E-01 

"MDAT,  NDAT"  PAIRS  -  1  1  1  2 
SPECIAL  COORDINATES  20.590  0.000 

PLOT  LIMITS  -  1160  2359  -.0236  -.0291 

HORIZONTAL  VERTICAL 

PLOT  SCALES  -  -.80000E-01  .40000E-01  .80000E-01  -.10000E+00  .40000E-01  .60000E-01 

"MDAT,  NDAT"  PAIRS  -  1  1  1  2 

SPECIAL  COORDINATES  20.590  0.000  0.000  0.000 

COORDINATES  AT  SHOT  EJECTION  -  -.20104E-01  -.66116E-01 

PLOT  LIMITS  -  1160  2359  -.0035  -.0090 

HORIZONTAL  VERTICAL 

PLOT  SCALES  -  -.20000E-01  .10000E-01  .30000E-01  -.20000E-01  .10000E-01  .30000E-01 

SUBROUTINE  CALCZ  -  SCALING  PARAMETERS 
2175  2342  3474 
1406  1490  2059 

COEFFICIENTS  OF  RECOIL  MOTION  CURVE 

-.337E-02  . 377E+01  .420E+00  -.157E+02  .235E+02  -.104E+02 

Z-ARRAY 

-.294E+01  -. 366E+02  -.486E+02  -.325E+01  -.366E+02  -.486E+02 
NZREF  =  2 

ZREF  =  0.00000 


"MDAT,  NDAT"  PAIRS  - 
WRITING  ARRAY  ON  FILE 

6 

2  2 

INDICES 

1  2 
ARE 

2 

6 

2 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

WRITING 

ARRAY 

ON 

FILE 

- 

INDICES 

ARE 

2 

1 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

WRITING 

ARRAY 

ON 

FILE 

- 

INDICES 

ARE 

2 

2 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

WRITING 

ARRAY 

ON 

FILE 

- 

INDICES 

ARE 

6 

1 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

WRITING 

ARRAY 

ON 

FILE 

- 

INDICES 

ARE 

6 

2 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

WRITING 

ARRAY 

ON 

FILE 

- 

INDICES 

ARE 

1 

1 

BR0512 

1 

1 

1 BRL 

TUBE 

HOR 

WRITING 

ARRAY 

ON 

FILE 

- 

INDICES 

ARE 

6 

1 

8R0512 

1 

1 

lBRI, 

TUBE 

HOR 

WRITING 

ARRAY 

ON 

FILE 

- 

INDICES 

ARE 

6 

2 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

WRITING 

ARRAY 

ON 

FILE 

- 

INDICES 

ARE 

1 

1 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

PLOTTING  IN  SUBROUTINE  UNITV 


"MDAT,  NDAT"  PAIRS  -  6  1 

SPECIAL  COORDINATES  20.590  0.000 

PLOT  LIMITS  -  1960  2159  -1.0000  .9900 

(0)MS(1) 

PLOT  SCALES  -  -.10000E+01  .26257E+72  .99000E+00  -.53397E-03  .26257E+72  -.11864E-03 


"MDAT,  NDAT"  PAIRS  -  6  2 

SPECIAL  COORDINATES  20.590  0.000 

PLOT  LIMITS  -  1960  2159  -1.0000  .9900 

(0)MS(1) 

PLOT  SCALES  -  -.10000E+01  .26257E+72  .99000E+00  -.31739E-02  .26257E+72  -.10133E-02 


"MDAT,  NDAT"  PAIRS  -  6  1  6  2 

SPECIAL  COORDINATES  20.590  0.000  0.000  0.000 

COORDINATES  AT  SHOT  EJECTION  -  -.11864E-03  -.24520E-02 

PLOT  LIMITS  -  1960  2159  -.0004  -.0002 

PLOT  SCALES  -  -.50000E-03  .50000E-03  .50000E-03  -.50000E-03 

WRITING  ARRAY  ON  FILE  -  INDICES  ARE  6  1  BR0512 

WRITING  ARRAY  ON  FILE  -  INDICES  ARE  6  2  BR0512 

WRITING  ARRAY  ON  FILE  -  INDICES  ARE  1  1  BR0512 


. 50000E-03  .50000E-03 

1  1  1BRL  TUBE  HOR 

1  1  1BRL  TUBE  HOR 

1  1  1BRL  TUBE  HOR 


FILTER  PARAMETERS 
NFILT  =  31  JTYPE 

EDGE  0.0000  .4000 
FX  1.00 
WTX  1.0 

FILTER  COEFFICIENTS 
1  .16166620E-04 

6  -. 18862205E-02 
11  . 19559239E-01 

16  0. 


2  NBANDS  = 


-.74131438E-04 

.32833366E-02 

-.29218184E-01 


.21102257E-03 

-.54084061E-02 

.44647104E-01 


-.49246632E-03 

.85345464E-02 

-.73733156E-01 


.  10096462E-02 
-. 13047 140E-01 
.  15615738E+00 


FILTER  PARAMETERS 

NFILT  =  127  FP  = 

.00500  TBWID  = 

.01541  DP  =  .01000 

DS  =  .01000 

FILTER  COEFFICIENTS 

1  -.57611818E-02 

-.  13120665E-02 

-. 14458869E-02 

-. 1577981  IE-02 

6  -. 18272702E-02 

-.  19407950E-02 

-. 20444951E-02 

-.21345667E-02 

11  -. 22685610E-02 

-. 23074309E-02 

-. 23235239E-02 

-. 23165943E-02 

16  -. 22201999E-02 

-. 21269670E-02 

-. 20019212E-02 

-.  18400413E-02 

21  -. 14040333E-02 

-. 11271255E-02 

-.80807552E-03 

-.45054632E-03 

26  .39177683E-03 

.88090117E-03 

. 14120564E-02 

•  19883328E-02 

31  .32498950E-02 

. 39438586E-02 

.46685324E-02 

•54340488E-02 

36  . 70536956E-02 

. 79036343E-02 

.87797060E-02 

.96719531E-02 

41  . 1 1502901E-01 

. 12431251E-01 

. 1336 1246E-01 

.  14291479E-01 

46  . 16127318E-01 

. 17023636E-01 

. 17900295E-01 

.  18750877E-01 

51  . 20359358E-01 

. 21 108628E-01 

.21814687E-01 

.  22474542E-01 

56  .23637514E-01 

.24135636E-01 

.  24573834E-01 

. 24947050E-01 

61  .25501145E-01 

. 25673888E-01 

.25780396E-01 

.  25814786E-01 

17054889E-02 
-.22110148E-02 
22818637E-02 
16423613E-02 
-.49316320E-04 
.  25922726E-02 
.62277021E-02 
.  10582451E-01 
.  15214561E-01 
.  19572419E-01 
.  23082659E-01 
•  25256526E-01 


o  PLOTTING  IN  SUBROUTINE  CALCZ 

r\> 

"MDAT,  MDAT"  PAIRS  -  2  1 

SPECIAL  COORDINATES  20.590  0.000 

PLOT  LIMITS  -  1160  2359  -9.0000  2.9900 

(0)MS(1)  (0)CM(1) 

PLOT  SCALES  -  -.90000E+01  .26257E+72  .29900E+01  -.41374E-01  .26257E+72  -.20370E-01 


"MDAT,  NDAT"  PAIRS  -  2  2 

SPECIAL  COORDINATES  20.590  0.000 

PLOT  LIMITS  -  1160  2359  -9.0000  2.9900 

(0)MS(1)  (0)CM(1) 

PLOT  SCALES  -  -.90000E+01  .26257E+72  .29900E+01  -.96551E-01  .26257E+72  .12158E-02 


"MDAT,  NDAT"  PAIRS  -  2  1  2  2 

SPECIAL  COORDINATES  20.590  0.000  0.000  0.000 

COORDINATES  AT  SHOT  EJECTION  -  -.20370E-01  -.72367E-01 

PLOT  LIMITS  -  1160  2359  -.0040  -.0095 

(0)CM(1)  (0)CM(1) 

PLOT  SCALES  -  -.20000E-01  .10000E-01  .20000E-01  -.20000E-01  .10000E-01  .20000E-01 


"MDAT,  NDAT"  PAIRS  -  3  1 

SPECIAL  COORDINATES  20.590  0.000 

PLOT  LIMITS  -  1160  2359  -9.0000  2.9900 

(0)MS(1)  (0)CM/MS( 1) 

PLOT  SCALES  -  -.10000E+02  .20000E+01  .40000E+01  -.15000E-01  .50000E-02  .15000E-01 

"MDAT,  NDAT"  PAIRS  -  3  2 

SPECIAL  COORDINATES  20.590  0.000 

PLOT  LIMITS  -  1160  2359  -9.0000  2.9900 

(0)MS(1)  (0)CM/MS( 1 ) 

PLOT  SCALES  -  -.90000E+01  .26257E+72  .29900E+01  -.56225E-01  .26257E+72  .39069E-01 

"MDAT,  NDAT"  PAIRS  -  3  1  3  2 

SPECIAL  COORDINATES  20.590  0.000  0.000  0.000 

COORDINATES  AT  SHOT  EJECTION  -  -.12763E-02  -.28196E-01 

PLOT  LIMITS  -  1960  2159  .0091  -.0032 

(o)cm/ms( 1 )  (o)cm/ms( 1) 

PLOT  SCALES  -  -.20000E-01  .10000E-01  .20000E-01  -.20000E-01  .10000E-01  .20000E-01 

g  ZREF  =  -7.52000 

"MDAT,  NDAT"  PAIRS  -  6  2  2  1  2  2 


WRITING 

ARRAY 

ON 

FILE 

-  INDICES 

ARE 

6 

2 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

WRITING 

ARRAY 

ON 

FILE 

-  INDICES 

ARE 

2 

1 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

WRITING 

ARRAY 

ON 

FILE 

-  INDICES 

ARE 

2 

2 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

WRITING 

ARRAY 

ON 

FILE 

-  INDICES 

ARE 

6 

1 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

WRITING 

ARRAY 

ON 

FILE 

-  INDICES 

ARE 

6 

2 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

WRITING 

ARRAY 

ON 

FILE 

-  INDICES 

ARE 

1 

1 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

WRITING 

ARRAY 

ON 

FILE 

-  INDICES 

ARE 

6 

1 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

WRITING 

ARRAY 

ON 

FILE 

-  INDICES 

ARE 

6 

2 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

WRITING 

ARRAY 

ON 

FILE 

-  INDICES 

ARE 

1 

1 

BR0512 

1 

1 

1BRL 

TUBE 

HOR 

PLOTTING  IN  SUBROUTINE  UNITV 


"MDAT,  NDAT"  PAIRS  -  6  1 

SPECIAL  COORDINATES  20.590  0.000 

PLOT  LIMITS  -  1960  2159  -1.0000  .9900 

(0)MS(1) 

PLOT  SCALES  -  -.10000E+01  .26257E+72  .99000E+00  -.53124E-03  .26257E+72  -.11726E-03 

"MDAT,  NDAT"  PAIRS  -  6  2 

SPECIAL  COORDINATES  20.590  0.000 

PLOT  LIMITS  -  1960  2159  -1.0000  .9900 

(0)MS(1) 

PLOT  SCALES  -  -.10000E+01  .26257E+72  .99000E+00  -.31721E-02  .26257E+72  -.10128E-02 


o 


"MDAT,  NDAT"  PAIRS  -  6  1  6  2 

SPECIAL  COORDINATES  20.590  0.000 

COORDINATES  AT  SHOT  EJECTION  -  -.11726E-03 

PLOT  LIMITS  -  1960  2159  -.0004  -.0002 

PLOT  SCALES  -  -.50000E-03  .50000E-03 

WRITING  ARRAY  ON  FILE  -  INDICES  ARE  6  1 

WRITING  ARRAY  ON  FILE  -  INDICES  ARE  6  2 

WRITING  ARRAY  ON  FILE  -  INDICES  ARE  1  1 


0.000  0.000 
-.24504E-02 

50000E-03  -.50000E-03 

BR0512 
BR0512 
BR0512 


.50000E-03  .50000E-03 

1  1  1BRL  TUBE  HOR 

1  1  1BRL  TUBE  HOR 

1  1  1 BRL  TUBE  HOR 


PLOTTING  IN  SUBROUTINE  CALCZ 

"MDAT,  NDAT"  PAIRS  -  2  1 

SPECIAL  COORDINATES  20.590  0.000 

PLOT  LIMITS  -  1160  2359  -9.0000  2  9900 

(0)MS(1)  (0)CM(1) 

PLOT  SCALES  -  -.90000E+01  .26257E+72  .29900E+01  -.35949F.-01  .26257E+72  -.19480E-01 


"MDAT,  NDAT"  PAIRS  -  2  2 

SPECIAL  COORDINATES  20.590  0.000 

PLOT  LIMITS  -  1160  2359  -9.0000  2.9900 

(0)MS(1)  (0)CM(1) 

PLOT  SCALES  -  -.90000E+01  .26257E+72  .29900E+01  -.72686E-01  .26257E+72  .20543E-02 


"MDAT,  NDAT"  PAIRS  -  2  1  2  2 

SPECIAL  COORDINATES  20.590  0.000  0.000  0.000 

COORDINATES  AT  SHOT  EJECTION  -  -.19480E-01  -.53931E-01 

PLOT  LIMITS  -  1160  2359  -.0023  -.0073 

(0)CM(1)  (0)CM(1) 

PLOT  SCALES  -  -.20000E-01  .10000E-01  .20000E-01  -.20000E-01  .10000E-01  .20000E-01 

"MDAT,  NDAT"  PAIRS  -  3  1 

SPECIAL  COORDINATES  20.590  0.000 

PLOT  LIMITS  -  1160  2359  -9.0000  2.9900 

(0)MS(1)  (0)CM/MS( 1 ) 

PLOT  SCALES  -  -.10000E+02  .20000E+01  .40000E+01  -.15000E-01  .50000E-02  .15000E-01 

"MDAT,  NDAT"  PAIRS  -  3  2 

SPECIAL  COORDINATES  20.590  0.000 

PLOT  LIMITS  -  1160  2359  -9.0000  2.9900 

(0)MS(1)  (0)CM/MS( 1 ) 

PLOT  SCALES  -  -.90000E+01  .26257E+72  .29900E+01  -.39433E-01  .26257E+72  .27812E-01 

g  "MDAT,  NDAT"  PAIRS  -  3  1  3  2 

SPECIAL  COORDINATES  20.590  0.000  0.000  0.000 

COORDINATES  AT  SHOT  EJECTION  -  -.19140E-02  -.22843E-01 

PLOT  LIMITS  -  1960  2159  .0063  -.0027 

(o)cm/ms( 1 )  (0)CM/MS( 1) 

PLOT  SCALES  -  -.20000E-01  .10000E-01  .20000E-01  -.20000E-01  .10000E-01  .20000E-01 
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D.2.6.  Sample  Plotted  Output 


These  plots  (Figures  D. 1 
Figures  D.l.  through  D.5,  the 
the  extrapolation  of  the  data 


through  D.25)  are  representative  of  the  ones  available  from  the  program, 
solid  line  refers  to  the  data  from  sensor  1  while  the  dashed  line  refers 
from  sensors  2  and  3  to  the  position  of  sensor  1. 
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Figure  D.l. 


Horizontal  Displacement , cm. 
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Figure  D.3,  Magnitude,  cm 


Figure  D,2.  Vertical  Displacement,  cm 


Figure  D.4,  Phase  Angle,  deg 


VERTICAL 

-O.IO  -0.06  -0.02  0.02  0.06 


Figure  D*5.  X-Y  Plot  Of  Displacement,  cm 
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Figure  D.6.  X-Y  Plot  Of  Sensor  1 
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Figure  D.7 


X-Y  Plot  Of  Sensor  1  Displacement,  Amplitude  Zeroed  At  Time  Of  Shot  Ejection 
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Figure  D.8. 


e'31  At  ZFEF 
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Figure  D.9.  e'^  ZREF 


Figure  D.ll. 


Horizontal  Displacement  At  ZREF 
=  .0,  cm 


Figure  D.12.  Vertical  Displacement  At  ZREF 
=  .  0  ,  cm 


Figure  D.13.  X-Y  Plot  Of  Displacement  At  ZREF  =  .0,  cm 
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Figure  D.14.  Horizontal  Velocity  At  ZREF  =  .0,  cm/ms  Figure  D.15.  Vertical  Velocity  At  ZFEF  =  .0,  cm/ms 


Figure  D.16.  JT-Y  Plot  Of  Velocity  At  ZREF  =  .0,  cm/ms 


ii.w.n  mi  tr«rtSE  IF  i.d  iBu  «7E 


O  o 


[ns] 


Figure  D.18.  e'32  At  ZREF  =  -  7.524  cm 
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Figure  D.19 


X-Y  Plot  Of  e' 


And  e' 


At  ZREF 


7.524  cm 
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Figure  D.20 
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Figure  D.21.  Vertical  Displacement  At  ZREF 
=  -  7.524  cm 


Figure  D.22.  J-Y  Plot  Of  Displacement  AT  ZREF 
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Figure  D.25 


X-Y  Plot  Of  e' 


And  e' 


At  ZREF 


7.524  cm 
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